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' Abstract. This is a survey of recent results on quantum ergodicity, specifically on the large 

, energy limits of matrix elements relative to eigenf unctions of the Laplacian. It is mainly 

^ ' devoted to QUE (quantum unique ergodicity) results, i.e. results on the possible existence 

Q , of a sparse subsequence of eigenfunctions with anomalous concentration. We cover the lower 

' bounds on entropies of quantum limit measures due to Anantharaman, Nonnenmacher, and 

^^.^ , Riviere on compact Riemannian manifolds with Anosov flow. These lower bounds give new 

' constraints on the possible quantum limits. We also cover the non-QUE result of Hassell in 

the case of the Bunimovich stadium. We include some discussion of Hecke eigenfunctions 

and recent results of Soundararajan completing Lindenstrauss' QUE result, in the context of 

'"j , matrix elements for Fourier integral operators. Finally, in answer to the potential question 

'^1^ ' 'why study matrix elements' it presents an application of the author to the geometry of 

r~| , nodal sets. 
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Quantum chaos on Riemannian manifolds (M, g) is concerned with the asymptotics of 
eigenvalues and orthonormal bases of eigenfunctions 

(1) Alpj = X'^jLfj, {^j,^k) = 5jk 

of the Laplacian A when the geodesic flow g^ : S*M S*M is ergodic or Anosov or in some 
other sense chaotic. Model examples include compact or flnite volume hyperbolic surfaces 
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Xr = r\H, where H is the hyperbohc plane and T C PSL{2,M) is a discrete subgroup. Of 
special interest are the arithmetic quotients when F is a discrete arithmetic subgroup. Other 
model example include Euclidean domains with ergodic billiards such as the Bunimovich 
stadium. The general question is, how does the dynamics of the geodesic (or billiard) flow 
make itself felt in the eigenvalues and eigenfunctions (or more general solutions of the wave 
equation), which by the correspondence principle of quantum mechanics must reproduce the 
classical limit as the eigenvalue tends to infinity? 

To avoid confusion, we emphasize that we denote eigenvalues of the Laplacian by They 
are usually viewed as energies Ej = Xj. Their square roots Xj are called the frequencies. 

The most basic quantities testing the asymptotics of eigenfunctions are the matrix elements 
{Aipj, ipj) of pseudo-differential operators relative to the orthonormal basis of eigenfunctions. 
As explained below, these matrix elements measure the expected value of the observable A in 
the energy state ipj. Much of the work in quantum ergodicity since the pioneering work of A. 
I. Schnirelman |Sh.l] is devoted to the study of the limits of {Aipj, ipj) as Xj oo. Although 
difficult to determine, these limits are still the most accessible aspects of eigenfunctions. A 
sequence of eigenfunctions is called ergodic or diffuse if the limit tends to /c.*» , crAdfiL where 

Og ill 

a A is the prinicipal symbol of A and dfiL is the Liouville measure. 

The first priority of this survey is to cover the important new results on the so-called 
scarring or QUE (quantum unique ergodicity) problem. Roughly speaking, the problem is 
whether every orthonormal basis becomes equidistributed in phase space with respect to 
Liouville measure as the eigenvalue tends to infinity, i.e. 



or whether there exists a sparse exceptional sequence which has a singular concentration. One 
series of positive results due to N. Anantharaman [S], N. Anantharaman-S. Nonnenmacher 
|ANj (see also [ANK] ) . and G. Riviere [Rivj fiows from the study of entropies of quantum 
limit measures and are based on a difficult and complex microlocal PDE analysis of the long 
time behavior of the wave group on manifolds with Anosov geodesic fiow. These results given 
lower bounds on entropies of quantum limit measures which imply that the limit measures 
must be to some extent diffuse. Another series of results due R. Soundararajan |Soundl] 
use L-function methods to complete the near QUE result of E. Lindenstrauss [LIND] for 
Hecke eigenfunctions in the arithmetic case. In the negative direction, A. Hassell |Has] uses 
microlocal methods to proved the long-standing conjecture that generic stadia are not QUE 
but rather have exceptional sequences of bouncing ball type modes which concentrate or scar 
along a Lagrangian submanifold. 

QUE is not the only interesting problem in quantum chaos. Another series of results in the 
arithmetic case are the remarkably sharp Luo-Sarnak asymptotics |LSt ILS2j of the variances 



for holomorphic Hecke eigenforms and their recent generalization to smooth Maass forms by 
Zhao |Zhj . Physicists (see e.g. |FP] ) have speculated that the variance is related to the clas- 
sical auto-correlation function of the geodesic flow. The Luo-Sarnak- Zhao results partially 
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confirm this conjecture in the special case of arithmetic surfaces and Hecke eigenfunctions, 
but correct it with an extra arithmetic factor. 

A final direction we survey is the applications of quantum ergodicity to problems in nodal 
geometry. Given the large amount of work on matrix elements {A(fj,(fj), it is natural for 
geometric or PDE analysts to ask how one can use such matrix elements to study classical 
problems on eigenfunctions such as growth and distributions of zeros and critical points. 
One of our themes is that for real analytic (M, g), logarithms log |v5j'(C)P of analytic contin- 
uations of ergodic eigenfunctions are asymptotically 'maximal pluri-subharmonic functions' 
in Grauert tubes and as a result their zero sets have a special limit distribution |Z5j . Thus, 
ergodicity causes maximal oscillation in the real and complex domain and gives rise to special 
distributions of zeros and (conjecturally) to critical points. 

In preparing this survey, we were only too aware of the large number of surveys and 
expository articles that already exist on the QUE and entropy problems, particularly |AN2t 
RNKl [CV3l ILIND21 [Sa?3] (and for earlier work [ZH [Z2l |M1 Ear] ) . To a large degree we closely 
follow the original sources, emphasizing intuition over the finer technical details. Although 
we go over much of the same material in what are now standard ways, we also go over 
some topics that do not seem as well known, and possibly could be improved, and we fill in 
some gaps in the literature. In particular, most expositions emphasize the quantum entropic 
uncertainty principle of [AN] due to its structural nature. So instead we emphasize the ideas 
of [X], which are less structural but in some ways more geometric. 

One of the novelties is the study in ^ of matrix elements {Fipj,(fj) of Fourier integral 
operators F relative to the eigenfunctions. Hecke operators are examples of such F but 
restrictions to hypersurfaces provide a different kind of F which have recently come up 
in quantum ergodic restriction theory |TZ3j . We give a rather simple result on quantum 
limits for Fourier integral operators which answers the question, "what invariance property 
do quantum limit measures of Hecke eigenfunctions possess?" . Hecke analysts have worked 
with a partial quasi-invariance principle due to Rudnick-Sarnak, and although the exact 
invariance principle may not simplify the known proofs of QUE it is of interest to know that 
one exists. We plan to put complete proofs in a future article. 

Another aspect of eigenfunctions is their possible concentration of eigenfunctions around 
closed hyperbolic geodesic. Mass concentration within thin or shrinking tubes around such 
geodesies has been studied by S.Nonnenmacher-A.Voros [NVj . Y. Colin de Verdiere-B. Parisse 
|CVPj . J. Toth and the author [TZ2] and outside of such tubes by N. Burq-M. Zworski |BZj 
and H. Christianson [Chr] . Very striking and surprising studies of analogous behavior in the 
'quantum cat map' setting are in [FNB[ [FN] . 

A few words on the limitations of this survey. Due to the author's lack of knowledge, the 
L-function or arithmetic methods in QUE will not be discussed in more than a cursory way. 
Rather we concentrate on the PDE or microlocal aspects of quantum ergodicity. By necessity, 
microfocal analytical proofs must make direct connections between spectrum and geodesic 
flow, and cannot by-pass this obstacle by taking a special arithmetic route. Moreover, we 
restrict the setting to Laplacians on Riemannian manifolds for simplicity of exposition, but 
it is just a special case of the semi-classical asymptotics (as the Planck constant tends to 
zero) of spectra of Schrodinger operators h'^A + V. Schrodinger eigenfunctions are equally 
relevant to physicists and mathematics. Another important setting is that of quantum 
'cat maps' on tori and other Kahler manifolds. In this context, the quantum dynamics 
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is defined by quantizations of symplectic maps of Kahler manifolds acting on spaces of 
holomorphic sections of powers of a positive line bundle. The mathematics is quite similar 
to the Riemannian setting, and results in the cat map setting often suggest analogues in 
the Riemannian setting. In particular, the scarring results of Faure-Nonnemacher- de Bievre 
|FNB] and the variance results of Kurlberg-Rudnick |KRj and Schubert |Schu3j are very 
relevant to the scarring and variance results we present in the Riemannian setting. 

Although we do not discuss it here, a good portion of the literature of quantum chaos is 
devoted to numerics and physicists' heuristics. Some classics are |Bert IHet IHe2j . The article 
of A. Barnett [Bar] (and in the cat map case |FNBl IFNj ) is a relatively recent discussion of 
the numerical results in quantum ergodicity. There is a wealth of phenomenology standing 
behind the rigorous results and heuristic proofs in this field. Much of it still lies far beyond 
the scope of the current mathematical techniques. 

Finally, we thank N. Anantharaman, A. Hassell, and S. Nonnenmacher for explanations 
of important points in their work. We have incorporated many of their clarifications in the 
survey. We also thank K. Burns, H. Hezari, J. Franks, E. Lindenstrauss, P. Sarnak, A. 
Wilkinson and J. Wunsch for further clarifications and corrections. Of course, responsibility 
for any errors that remain is the author's. 

1. Wave equation and geodesic flow 

Classical mechanics is the study of Hamiltonians (real valued functions) on phase space 
(a symplectic manifold) and their Hamiltonian flows. Quantum mechanics is the study of 
Hermitian operators on Hilbert spaces and the unitary groups they generate. The model 
quantum Hamiltonians we will discuss are Laplacians A on compact Riemannian manifolds 
{M,g) (with or without boundary). Throughout we denote the dimension of M by d = 
dimM. 

1.1. Geodesic flow and classical mechanics. The classical phase space in this setting is 
the cotangent bundle T*M of M, equipped with its canonical symplectic form J2i A rf^j. 

The metric defines the Hamiltonian H{x,C,) = \C,\g = \JY1%=i 9^K^)iiij T*M, where 
~ af")' t^'*"'] inverse matrix to [gij]. We denote the volume density of {M,g) 

by dVol and the corresponding inner product on L'^{M) by (/, g). The unit (co-) ball bundle 
is denoted B*M = {(x,0 : I'^l < !}• 

The classical evolution (dynamics) is given by the geodesic flow of {M,g), i.e. the Hamil- 
tonian flow ^f* of H on T*M. By definition, g^{x,^) = {xt,C,t), where {xt,^t) is the terminal 
tangent vector at time t of the unit speed geodesic starting at x in the direction ^. Here and 
below, we often identify T*M with the tangent bundle TM using the metric to simplify the 
geometric description. The geodesic flow preserves the energy surfaces {H = E} which are 
the co-sphere bundles S*^M. Due to the homogeneity of H, the flow on any energy surface 
{H = E} is equivalent to that on the co-sphere bundle S*M = {H =1}. 

We define the Liouville measure fiL on S* M to be the surface measure dfi^ = ^f^- induced 
by the Hamiltonian H = \^\g and by the symplectic volume measure dxdC, on T*M. The 
geodesic flow (like any Hamiltonian flow) preserves dxdC,. It also preserves Liouville measue 
on S*M. 
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1.2. Ergodic and Anosov geodesic flows. The geodesic flow is called ergodic if the 
unitary operator 

(2) V'a{p) = a{g\p)) 

on L'^{S*M, diiL) has no invariant functions besides constant functions. Equivalently, any 
invariant set E d S*M has either zero measure or full measure. 

A geodesic flow (7* is called Anosov on on S*M if the tangent bundle TS*M splits into g* 
invariant sub-bundles E^{p) © E^{p) © M.Xh{p) where E"^ is the unstable subspace and E'^ 
the stable subspace. They are deflned by 

\\dg^v\\ < Ce-^^\\v\\, 'rfveE',t>0, 

\\dg^v\\ < Ce^^\\v\\, WveE'',t<0. 

The sub-bundles are integrable and give stable, resp. unstable foliations W^jW^. The 
leaves through x are denoted by W^i^x), Thus, the geodesic flow contracts everything 

exponentially fast along the stable leaves and expands everything exponentially fast along 
the unstable leaves. Anosov geodesic flows are ergodic. We refer to pCl] for background. 
The unstable Jacobian J"(p) at p is deflned by 

(3) J"(p) = det^-l(^,^)). 
Define 

A = -suppgs.^jlog||J"(p)|| 

(4) 

= inf^eA^, \ J g,^.j log J''{p)dp{p)\ > 0. 

If {M,g) is of dimension d and constant curvature —1 then A = d — 1. Here, Aij is the set 
of invariant probability measures for g^. 

The maximal expansion rate of the geodesic flow is deflned by 

(5) Amax=lim^log sup \\dgp\\. 

The derivative dg^ may be expressed in terms of Jacobi flelds (see [Klj . Lemma 3.1.17). 
In the hyperbolic (Anosov case), the stable/unstable sub-bundles are spanned by unstable 
Jacobi flelds Y^, . . . , Y]^_-^ whose norms grow exponentially as t — > +00 and decay exponen- 
tially a.st —>■ —00, resp. stable Jacobi flelds Yi, . . . , FjL^ whose norms behave in the opposite 
fashion. The unstable Jacobian is the norm \ \Y^ A ■ ■ ■ A 

1.3. Eigenfunctions and eigenvalues of A. The quantization of the Hamiltonian H is 
the square root \/A of the positive Laplacian, 

^ L ^ «i ^ 



V9 t]^i dxi dxj 

of {M,g). Here, g = det[gij]. The eigenvalue problem on a compact Riemannian manifold 



k 



is dual under the Fourier transform to the wave equation. Here, {^pj} is a choice of orthonor- 
mal basis of eigenfunctions, which is not unique if the eigenvalues have multiplicities > 1. 
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The individual eigenfunctions are difficult to study directly, and so one generally forms the 
spectral projections kernel, 

(6) E{\,x,y)= ipj{x)ipj{y). 

Semi-classical asymptotics is the study of the A oo limit of the spectral data {ipj, Xj} or 
of E{\x,y). 

1.4. Eigenvalues and Planck's constant. Although it is only a notational issue in this 
setting, it is conceptually useful to identify inverse frequencies with Planck's constant, i.e. 
to define 

(7) n, = xj\ 

The eigenvalue problem then takes the form of a Schrodinger equation 

h'^jAifj = (fj. 

Thus the high frequency limit is a special case of the semi-classical limit h 0. 

1.5. Eigenvalue Multiplicities. The multiplicity of an eigenvalue is the dimension 

(8) m{Xj) = dimker(A - A^^) 

of the eigenspace. As will be discussed in ^ the general Weyl law gives the bound m(Aj) < 
CA^~^ in dimension d. In the case of negatively curved manifolds, the bound has been 
improved by a (logAj)"^ factor. Thus it is possible for sequences of eigenvalues to have 
the huge multiplicities A^~^(log Aj)^^. When such multiplicities occur, it is possible to build 
up superpositions of eigenfunctions which have special concentration properties (see 
Hassell's scarring theorem is based on the non-existence of such multiplicities for certain 
sequences of eigenvalues (§7]). 

In fact such multiplicities do occur in the quantum cap map setting for sparse sequences 
of eigenvalues and are responsible for the very strange behavior of eigenstates in that case 
|FNBt IFN] . It is unknown if such sequences exist in the Riemannian Anosov setting. This 
is one of the fundamental obstacles to understanding whether or how scarring occurs in the 
Riemannian case. 

1.6. Quantum evolution. Quantum evolution is given by the wave group 

which in a rigorous sense is the quantization of the geodesic flow $*. It is generated by 
the pseudo-differential operator V^, as defined by the spectral theorem (it has the same 
eigenfunctions as A and the eigenvalues Xj). The (Schwartz) kernel of the wave group on a 
compact Riemannian manifold can be represented in terms of the spectral data by 

(9) U\x,y) = Y,e''^'^A^WAy)^ 

j 

or equivalently as the Fourier transform e**^(ii?(A, x, y) of the spectral projections. Hence 
spectral asymptotics is often studied through the large time behavior of the wave group. 
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Evolution can be studied on the level of 'points' or on the level of observables. Evolution of 
points in classical mechanics gives the orbits or trajectories of the geodesic flow (i.e. the pa- 
rameterized geodesies). Evolution ipt = U^ip of wave functions gives the 'Schrodinger picture' 
of quantum mechanics. Eigenf unctions arise in quantum mechanics as stationary states, i.e. 
states ip for which the probability measure \ip(t,x)\'^dvol is constant where ipit^x) = U^tp{x) 
is the evolving state. This follows from the fact that 

(10) U'cpk = e^*^Vfc 

and that |e**^'- | = 1. By unitarity f/** = f/"*. 

As an alternative one could define the quantum evolution as in jXN] to be the semi-classical 
Schrodinger propagator 

(11) Ui := e"'*'^^. 

This replaces the homogeneous Fourier integral operator by a semi-classical Fourier inte- 
gral operator. The theories are quite parallel; we mainly keep to the homogeneous theory in 
this survey (See |DSj| and |EZj ) for background on these notions). 



1.7. Observables. In the classical setting, observables are (real- valued) functions on T*M 
which we usually take to be homogeneous of degree zero (or as zeroth order symbols). In 
the quantum setting, observables are ifjDO^s (pseudo differential operators) of all orders; we 
often restrict to the subalgebra \Ef° of ipDO^s of order zero. We denote by \E'"^(M) the 
subspace of pseudodifferential operators of order m. The algebra is defined by constructing 
a quantization Op from an algebra of symbols a G S"^{T*M) of order m (polyhomogeneous 
functions on T*M\0) to \I''". A function on T*M\0 is called polyhomogeneous or a classical 
symbol of order m, a G S^f^g, if it admits an asymptotic expansion, 

oo 

0, where ak{x,T^) = r'^ttkix,^), (|^| > 1) 

j=0 

for some m (called its order). The asympotics hold in the sense that 

R 



r- cm-R-1 

j=0 



where a G S''^ if sup(l + \C,\y~''\D^D'^a{x,^)\ < +oo for all compact set and for all a,f],j. 
There is a semi-classical analogue in the h setting where the complete symbol afi{x, ^ ) depends 
on ^ in a similar poly-homogeneous way. 



R 

an ~ h^"^ hPttr, 



j=0 

We refer to |DSj[ IGStt lEZt IDuij for background on microlocal analysis (or semi-classical 



analysis), i.e. for the theory of pseudo-differential and Fourier integral operators. 

The main idea is that such observables have good classical limits. The high frequency 
limit is mirrored in the behavior of the symbol as |^| oo. 

An useful type of observable is a smooth cutoff Xu{x,^) to some open set U C S*M (it 
is homogeneous so it cuts off to the cone through U in T*M). That is, a function which is 
one in a large ball in U and zero outside U . Then Op{xu) is called a microlocal cutoff to 
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U. Its expectation values in the states give the phase space mass of the state in U or 
equivalently the probabihty amphtude that the particle is in U . (The modulus square of the 
amplitude is the probability). 

There are several standard useful notations for quantizations of symbols a to operators A 
that we use in this survey. We use the notation Op{a) or a interchangeably. 

1.8. Heisenberg uncertainty prinicple. The Heisenberg uncertainty principle is the 
heuristic principle that one cannot measure things in regions of phase space where the product 
of the widths in configuration and momentum directions is < h. 

It is useful to microlocalize to sets which shrink as h ^ using h (or A)-dependent 
cutoffs such as Xu{^^{x,^)). The Heisenberg uncertainty principle is manifested in pseudo- 
differential calculus in the difficulty (or impossibility) of defining pseudo-differential cut-off 
operators Op{xu{^^{x, ^)) when S > ^. That is, such small scale cutoffs do not obey the usual 
rules of semi-classical analysis (behavior of symbols under composition). The uncertainty 
principle allows one to study eigenfunctions by microfocal methods in configuration space 
balls B{xo,Xj^) of radius h (since there is no constraint on the height) or in a phase space 

ball of radius XJ^^^- 

1.9. Matrix elements and Wigner distributions. These lectures are mainly concerned 
with the matrix elements {Op{a)ipj, ipk) of an observable relative to the eigenfunctions. The 
diagonal matrix elements 

(12) p^{Op{a)):^{Op{a)^,,^,) 

are interpreted in quantum mechanics as the expected value of the observable Op{a) in the 
energy state ipj. The off-diagonal matrix elements are interpreted as transition amplitudes. 
Here, and below, an amplitude is a complex number whose modulus square is a probability. 

If we fix the quantization a — > Op{a), then the matrix elements can be represented by 
Wigner distributions. In the diagonal case, we define Wk G V'{S*M) by 

(13) / adWk :^ {Op{a)(fk,(fk)■ 
Js*M 

1.10. Evolution of Observables. Evolution of observables is particularly relevant in quan- 
tum chaos and is known in physics as the 'Heisenberg picture'. 

The evolution of observables in the Heisenberg picture is defined by 

(14) at{A) := U^AU'^ A e *'"(M). 

and since Dirac's Principles of Quantum Mechanics, it was known to correspond to the 
classical evolution 

(15) V\a) := a o 

of observables a G C°°{S*M). Egorov's theorem is the rigorous version of this correspon- 
dence: it says that at defines an order-preserving automorphism of ^*(M), i.e. at{A) e 
^'"(M) ff ^ e ^™(M), and that 

(16) c7u,AU*{x,0 = ^a{9\x,0) := V^*(^a), {x,0 ^ T*M\0. 

This formula is almost universally taken to be the definition of quantization of a fiow or map 
in the physics literature. 
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1.11. Hadamard parametrix. The link between spectral theory and geometry, and the 
source of Egorov's theorem for the wave group, is the construction of a parametrix (or WKB 
formula) for the wave kernel. For small times t, the simplest is the Hadamard parametrix. 



oo 



(17) U\x,y)^ e^'^^'^^'y^-''^Y.^u{x,y)e^~^de (t<inj(M,^)) 

where r{x,y) is the distance between points, Uo{x,y) = Q^^{x,y) is the volume 1/2-density, 
inj{M,g) is the injectivity radius, and the higher Hadamard coefficients are obtained by 
solving transport equations along geodesies. The parametrix is asymptotic to the wave 
kernel in the sense of smoothness, i.e. the difference of the two sides of (fT7|) is smooth. The 
relation (fTTl) may be iterated using (?7*)™ = t/*"* to obtain a parametrix for long times. For 
manifolds without conjugate points, the Hadamard parametrix is globally well defined on 
the universal cover and can then be summed over the deck transformation group to obtain 
a parametrix on the quotient. 

Egorov's theorem (ITB]) may be proved by explicitly writing out the integral representation 
for U^^AU^ using 0171) and by applying the stationary phase method to simplify the integral. 

1.12. Ehrenfest time. The aim in quantum chaos is to obtain information about the high 
energy asymptotics as A-, ^ oo of eigenvalues and eigenfunctions by connecting information 
about f/* and g^. The connection comes from comparing ([9]) and f|T7|) . or by using Egorov's 
theorem f|T6l) . But to use the hypothesis that g^ is ergodic or chaotic, one needs to exploit the 
connection as — >■ and t —>■ oo. The difficulty in quantum chaos is that the approximation 
of by g^ is only a good one for t less than the Eherenfest time 

(18) T, = i2i^. 

where Amax is defined in ([5]). 

Roughly speaking, the idea is that the evolution of a well constructed "coherent" quantum 
state or particle is a moving lump that "tracks along" the trajectory of a classical particle up 
to time Th and then slowly falls apart and stops acting like a classical particle. Numerical 
studies of long time dynamics of wave packets are given in works of E. J. Heller [Het IHe2] 
and rigorous treatments are in Bouzouina- Robert |BouRj . Combescure- Robert |CRj and 
Schurbert |Schu3] . 

The basic result expressed in semi-classical notation is that there exists F > such that 

(19) \\U^'Op{a)Ui-Op{aog')\\<Che'^. 

Such an estimate has long been known if F is not specified in terms of the geodesies. It is 
implicit in |Be] and used explicitly in [Z4], based on long time wave equation estimates of 
Volovoy. The best constant should be F = Amax- 

Thus, one only expects good joint asymptotics as ^ — 0, T ^ cx3 for t < Te- As a result, 
one can only exploit the approximation of by g^ for the relatively short time Te- 

1.13. Modes, quasi-modes and coherent states. This survey is mainly concerned with 
eigenfunctions. But one of the main tools for studying eigenfunctions is the construction 
of quasi-modes or approximate eigenfunctions. Lagrangian and coherent states are closely 
related to quasi-modes. We refer to [Duit ICV2t \B.B\ IRaj IRa2j for background. 
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A quasi-mode {ipk} of order zero is a sequence of L^-normalized functions satisfying 

(20) ||(A-^,)V'fc||L2 = 0(1), 

for a sequence of quasi-eigenvalues fik- By the spectral theorem, there must exist true 
eigenvalues in the interval [fik — 6, fi^ + 6] for some 6 > 0. Moreover, if E^^s denotes the 
spectral projection for the Laplacian corresponding to this interval, then 

\\Ek,5^k - ^k\\L2 = Oik"^). 
A more general definition is that a quasi-mode of order s is a sequence satisfying 

(21) ||(A-/i,)^fc|U2 = 0(/in- 
Then 

This definition allows for very weak versions of quasi- modes, e.g. quasi-modes of order — 1. 
They carry no information about eigenvalues since there always exist eigenvalues in intervals 
of width one. However, they do carry some information on eigenf unctions. 

In references such as |Duij . quasi-modes are constructed in the form of oscillatory integrals 
or Lagrangian states, 

^^(x) = [ e'^^(^'«)a(x,Of^^- 

The phase generates the Lagrangian submanifold 

As = {{x,d^S) : d^S{x,^) = 0} C T*M. 

If A C S'*M is a closed Lagrangian submanifold invariant under the geodesic flow then ipx 
is a quasi-mode of order —1. In favorable circumstances, one can further determine the 
amplitude so that ||(— A — A^)'?/'a|| = 0{\~'^) and construct a high order quasi-mode. Such 
a quasi-mode is the 'quantization' of A and the associated sequence concentrates on A. 

A coherent state is a special kind of Lagrangian state, or more precisely isotropic state, 
which quantizes a 'point' in phase space rather than a Lagrangian submanifold. The sim- 

plest examples are the states (p^°'^"{x) = /i^"/^es^^'^o^e on M", which are localized 

to the smallest possible volume in phase space, namely a unit cube of volume h around 
(xo,^o)- They are Lagrangian states with complex phases associated to positive Lagrangian 
submanifolds of complexified phase space. 

1.14. Riemannian random waves. Eigenf unctions in the chaotic case are never similar 
to such quasi-modes or Lagrangian states. Rather, they are conjectured to resemble random 
waves. This heuristic was first proposed by M. V. Berry [Berj . who had in mind random 
Euclidean plane waves. The specific model discussed here was used by the author to show 
that eigenfunctions in the ergodic case behave to some degree like random orthonormal bases 
of spherical harmonics. We use the term Riemannian random waves for this model. 

We decompose the spectrum of a compact Riemannian manifold (M, g) into intervals 
Jat = [N, A^ -|- 1] and define the finite dimensional Hilbert spaces Tii^ as the space of linear 
combinations of eigenfunctions with frequencies in an interval In := [A^, A^ -|- 1]. The precise 
decomposition of M into intervals is not canonical and any choice of intervals of finite width 
would work as long as the set of closed geodesies of (M, g) has measure zero. We will not 
discuss the other case since it has no bearing on quantum ergodicity. We denote by {v?Arj}^=i 
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an orthonormal basis of Tii^^ where = dimT-^/^^. We endow the real vector space Tii^ 
with the Gaussian probabihty measure 7Ar defined by 

(22) ^^(s)=(^\ e-'^-I^I'dc, ^P = J2cjVNj, dN = dimni^. 

V ^ / j=i 

Here, dc is (iAr-dimensional real Lebesgue measure. The normalization is chosen so that 
E^j^ {'ip^Tp) = 1, where E^^^ is the expected value with respect to •jn- Equivalently, the dN 
real variables Cj (j = l,...,diy) are independent identically distributed (i.i.d.) random 
variables with mean and variance ^j^; i.e., 

UN 

We note that the Gaussian ensemble is equivalent to picking ipN ^ 'Hi^ at random from the 
unit sphere in Hj^ with respect to the inner product. 

Numerical results confirming some features of the random wave model are given in [HRj . 

1.15. Hyperbolic plane. It is useful to have one concrete example in mind, and we will use 
the hyperbolic plane 7i or disc D as a running model. In the disc model D = {z G C, |z| < 1}, 
the hyperbolic metric has the form ds'^ = j-^^^j^. Its isometry group is G = PSU{1, 1); the 

stabilizer of is ~ 5*0(2) and thus D ~ G/K. Without comment we will also identify D 
with the upper half plane and G with PSL{2,M.). In hyperbolic polar coordinates centered 
at the origin 0, the (positive) Laplacian is the operator 

,9 1 92 
—A = — — ■ + cothr— — h 



dr'^ dr sinh^ r dO"^ 

The distance on D induced by the Riemannian metric will be denoted dj^. We denote the 
volume form by d{z). 

The unit tangent bundle SY) of the hyperbolic disc D may be identifed with the unit 
cosphere bundle 5**0 by means of the metric. We further identify 5*0 = PSU{1, 1), since 
PSU{1, 1) acts freely and transitively on 5*0. We identify a unit tangent vector {z,v) with 
a group element g if g ■ (0, (1, 0)) = {z, v). 

We denote hj B = {z & C,\z\ = 1} the boundary at infinity of D. Then we can also 
identify 5*0 = D x 5. Here, we identify {z,b) eD x B with the unit tangent vector {z,v), 
where v G S'zD is the vector tangent to the unique geodesic through z ending at b. 

The geodesic fiow g^ on SD is defined by g^{z,v) = (7^(t), 7^(t)) where 'jvit) is the unit 
speed geodesic with initial value {z,v). The space of geodesies is the quotient of SD by 
the action of (7*. Each geodesic has a forward endpoint b and a backward endpoint b' in B, 
hence the space of geodesies of D may be identified with B x B \ A, where A denotes the 
diagonal in _B x B: To (6', b) E B x B \ A there corresponds a unique geodesic ■jb/^f, whose 
forward endpoint at infinity equals b and whose backward endpoint equals b'. We then have 
the identification 

S'D = {BxB\A)xm. 

The choice of time parameter is defined - for instance - as follows: The point (6', b, 0) is by 
definition the closest point to on 7;,/.;, and {b',b,t) denotes the point t units from {b,b',0) 
in signed distance towards b. 
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1.15.1. Dynamics and group theory ofG = SL{2, M). The generators of sl{2, M) are denoted 
by 

"=(::)■-(::)■-(;: 

We denote the associated one parameter subgroups by A,A_,K. We denote the rais- 
ing/lowering operators for i^-weights by 

(23) E+ = H + iV, E- = H- iV. 

The Casimir operator is then given by 4 = + V^ — W"^; on invariant functions, the 
Casimir operator acts as the laplacian A. We also put 

and denote the associated subgroups by N,N_. In the identification 5*0 = PSL{2,'R) the 
geodesic flow is given by the right action g^{g) = gat of A, resp. the horocycle flow (/i")ugk 
is defined by h^{g) = gUu, where 
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The closed orbits of the geodesic flow g^ on V\G are denoted {7} and are in one-to-one 
correspondence with the conjugacy classes of hyperbolic elements of V. We denote by G^, 
respectively F^, the centralizer of 7 in G, respectively F. The group F^ is generated by an 
element 70 which is called a primitive hyperbolic geodesic. The length of 7 is denoted > 
and means that 7 is conjugate, in G, to 

(24) a. 



•1 







If 7 = 7q where 70 is primitive, then we call L^^ the primitive length of the closed geodesic 
7- 

1.15.2. Non-Euclidean Fourier analysis. Following [Hel], we denote by (-2,6) the signed 
distance to of the horocycle through the points z G D, 6 G -B. Equivalently, 

where -Pd(-2, &) is the Poisson kernel of the unit disc. We denote Lebesgue measure on B by 
\dh\. We then introduce the hyperbolic plane waves e*^2+*^)<^'*), the analogues of the Euclidean 
plane waves e*^^''^^ 

The non- Euclidean Fourier transform T : GcD ^ C{B x M_,_) is defined by 

J^u{r,b)= [ e^'^-'''^^''^^u{z)dVol{z). 
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The inverse Fourier transform is given by 

jr-i^(^) = f f e^-2+i'''>^^^b^g{r,b)rtanh{2nr)dr\db\. 



The integration measure is the Plancherel measure. extends to an isometry JF : L^(D, dV) — 
L'^{B X R+,rtanh(27rr)rfrc?6). 

1.15.3. Wave kernel on D and wave kernel on Xp. The wave group [/* = e**^ of the positive 
Laplacian on D is thus given in terms of Fourier analysis by (distribution) integral 



U\z,w) = / / e'(3+-')*e(5+-)(-.^)e(-5+-)Kf')rtanh(27rr)dr|d6|. 



The integral over B may be eliminated by the stationary phase method to produce an integral 
over M alone; it is the Hadamard parametrix (fT7|) in this case. 

Now let r C G be a discrete subgroup. In applications to quantum chaos we assume the 
quotient Xr is compact or of finite 

The wave kernel on the quotient Xr is obtained by automorphizing the wave kernel on D: 

(25) U'^{z,w) = Y,U'Az,iw). 

7er 

This representation is one of the inputs into the Selberg trace formula. The elements of F 
are grouped into conjugacy classes [7] G [F]. The conjugacy classes correspond to closed 
geodesies of Xp, i.e. periodic orbits of (?* on S'*Xr. 
Let 

(26) er(T) = #{[7]:L,<T}. 
The prime geodesic theorem asserts that 



(27) er(T) ^ , 

where htop is the topological entropy of g*. In fact, htop = 1 for the hyperbolic case. The 
exponential growth of the length spectrum reflects the exponential growth of the geodesic 
flow. 

The Ehrenfest time (1181) is implicit in the exponential growth rate of 0r(T). 

1.15.4. Representation theory of G and spectral theory of A. Let F C G be a co-compact 
discrete subgroup, and let us consider the automorphic eigenvalue problem on G/K: 

( A^ = (i+r2)^, 

(28) <^ 

I ^p{'jz) = (f{z) for all 7 G F and for all z. 

The solutions are the eigenfunctions of the Laplacian on the compact surface Xp = r\ G /K. 
A standard notation for the eigenvalues is = s(l — s) where s = | + ir. 

Eigenfunctions of the Laplacian are closely connected with the representation theory of G 
on G/r. We briefly describe the representation theory since the problems in quantum chaos 
have analogues in the discrete series as well as the unitary prinicipal series. 
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In the compact case, we have the decomposition into irreducibles, 

5 oo oo oo 

L2(r\G') = 0C,,.e07',,.e /^rMP+e //r(m)p-, 

j=l j=0 m=2, m even m=2,m even 

where Cir^ denotes the complementary series representation, respectively Vir^ denotes the 
unitary principal series representation, in which the Casimir operator —Q = —{H'^+V'^ — W'^) 
equals Sj(l — Sj) = \ + r^- In the complementary series case, irj G M while in the principal 
series case irj G iW^. The irreducibles are indexed by their i^- invariant vectors {(fiirj}, which 
is assumed to be the given orthonormal basis of A-eigenfunctions. Thus, the multiplicity of 
Vir is the same as the multiplicity of the corresponding eigenvalue of A. 

Further, denotes the holomorphic (respectively anti-holomorphic) discrete series rep- 
resentation with lowest (respectively highest) weight m, and jivijn) denotes its multiplicity; 
it depends only on the genus of Xp. We denote by (j = 1, • • • , l^^vijn)) a choice of or- 
thonormal basis of the lowest weight vectors of iXT{m)V^ and write iJ,r{m)T>:^ = (Bj=^^^T>^ j 
accordingly. 

There is a a direct integral decomposition for co-finite subgroups such as F = PSL{2, Z) 
or congruence subgroups. The non-compactness gives rise to a continuous spectral subspace 
of Eisenstein series and a discrete spectral subspace of cuspidal eigenfunctions (which is only 
known to be non-trivial in the case of arithmetic F). 

1.15.5. Helgason Poisson formula. The Fourier transform of an function on D is an 
fimction on 5 x ]R_|_. There is an extension of the FoTiricr transform and the inversion formula 
to tempered distributions. We only consider the case of F-automorphic eigenfunctions where 
F is co-compact. One then has J-^Pivj — dTirj{b) <S) Srj{r), where dTir^ is a distribution on the 
ideal boundary B. The inversion formula is Helgason's Poisson formula, 

(29) ^^r{z) = f e(i+-)<^'''>r,,,^,,(d6), 

Jb 

for all z G D. The kernel e*^^^*''^^^'''^ = P^^^^\z,b) is called the generalized Poisson kernel. 
The distribution 

(30) T,r,^Jdb) = J2(^nmdb\. 

neZ 

is called the boundary value of (fir and is obtained from the Fourier expansion 

(31) Vir{z) = y^an(r)$r-,n(^), 

of (fir in the disc model in terms of the generalized spherical functions, 

(32) S+ir){z,b) ^ J2 $,,„(^)6", b G B. 

nGZ 

Equivalently, the $r,n are the joint eigenfunctions of A and of K. 
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When ipir- is a F-invariant eigenfunction, the boundary values Tir^ (db) have the following 
invariance property: 

(33) 

This follows by the identities 

(34) {g-z,g-b) = {z,b) + {g-0,g-b), 
which implies 

(35) P^{gz,gb)\d{gb)\=Pj,{z,b)\db\. 

An interesting heuristic related to Berry's random wave hypothesis (see §1.14p is the 
conjecture that the Fourier coefficients a„(r) in fl5P]) should behave like independent Gaussian 
random variables of mean zero and variance one. This conjecture was stated explicitly and 
tested numerically by Aurich-Steiner in |AS] . Related tests of the random wave model for 
hyperbolic quotients are in Hejhal-Rackner [HRj . 

Otal [D] and Schmid [Schj have shown that Tir.{db) is the derivative of a Holder C2 
continuous function Fi^. on B. Since its zeroth Fourier coefficient is non-zero, Ti^^{db) is not 
literally the derivative of a periodic function, but it is the derivative of a function Fi^.^ on M 
satisfying Fir^{9 + 2n) = Fir^{9) + Cj for all 9 e R. Recall that for < 5 < 1 we say that 
a 27r-periodic function F : M ^ C is 5-H51der if \F{9) - F{9')\ < C\9 - 9'\^. The smallest 
constant is denoted and A5 denotes the Banach space of 5-Holder functions, up to 

additive constants. 

1.15.6. Boundary values and representation theory |Z3j . The distributions dTi^. have a nat- 
ural interpretation in representation theory. We define Cir^ G P'(F\P5'L(2, R)) such that 

(36) e(5+-.)(-''')T,,^ [db) = Cir, {z, b)P{z, b)db. 

The distribution Cir- is horocyclic- invariant and F-invariant. It may be expanded in a K- 
Fourier series, 

and it is easily seen (cf. ^Z2J) that ,0 = Vir- and that v^ir ,n is obtained by apply- 
ing the nth normalized raising or lowering operator (Maass operator) = H ± iV to 

(fir-. More precisely, one applies {E^)^ and multiplies by the normalizing factor (32ir ,n = 
^1 



{2irj+l±2n)-(2irj+l±2) ' 



2. WEYL law AND LOCAL WEYL LAW 



We now return to the general case. A basic result in semi-classical asymptotics is Weyl's 
law on counting eigenvalues: 

(37) N{\) = j^{j:\,<\} = ^Vol{M,g)\' + R{X), where R{\) = O {X"-^) . 
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Here, \Bd\ is the Euclidean volume of the unit ball and Vol{M,g) is the volume of M with 
respect to the metric g. Weyl's law says that 

voim, < A) 



(38) TrE, 



(27r)' 



where Vol is the symplectic volume measure relative to the natural symplectic form X]j=i dxjA 

d^j on T*M. Thus, the dimension of the space where H = ^/A is < A is asymptotically the 
volume where its symbol < A. 

The growth of the remainder term depends on the long time behavior of (7*. It is sharp 
on the standard sphere, where all geodesies are periodic. By a classical theorem due to 
Duistermaat-Guillemin |D.G] (in the boundaryless case) and to Ivrii (in the boundary case), 

-R(A) = o(A°'~^), when the set of periodic geodesies has Liouville measure zero. 

The remainder is then of small order than the derivative of the principal term, and one has 
asymptotics in shorter intervals: 

(39) N{[X, A + 1]) = #{j : A, e [A, A + 1]} = n^Vo/(M, g)X'-' + o{X'-'). 

Then the mean spacing between the eigenvalues in this interval is ~ CdV ol{M, g)^^ \^^'^^^\ 
where Cd is a constant depending on the dimension. 

In the case of compact Riemannian manifolds of negative curvature, a sharper estimate of 
the remainder is possible: 

(40) «(^) = 0(^)- 

This remainder is proved using fl25l) and the exponential growth rate in fl27|) . The estimate 
fHOj) was proved by Selberg in the case of compact hyperbolic quotients and was generalized 
to all compact Riemannian manifolds without conjugate points by Berard [Bej. The loga- 
rithm in the remainder is a direct outcome of the fact that one can only use the geodesic 
approximation to f/* up to the Ehrenfest time Te ffTSl) . 

This estimate has not been improved in fifty years, and there are no better results in 
the constant curvature case than in the general negatively curved case. The remainder 
estimate does not rule out the implausible-seeming scenario that in the interval [A, A + 1] 
there are only log A distinct eigenvalues with multiplicities In fact, such implausibly 

large multiplicities do occur for a sparse set of 'Planck constants' in the analogous case of 
quantizations of hyperbolic toral automorphisms (see [FNB] ). The issue of possibly high 
multiplicity or clustering of eigenvalues is important in Hassell's scarring result as well as 
the Anantharaman-Nonnenmacher results on entropies. 

An important generalization is the local Weyl law concerning the traces TrAE{X) where 
A e ^''^(M). It asserts that 

(41) = ^^"^^dA A'^ + O(A'^-i). 



There is also a pointwise local Weyl law: 
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(42) 5^|y.,(x)r = -i-^|i?^|A'^ + i?(A,x), 

A,<A ^ ^ 

where R{X, x) = 0{X'^~^) uniformly in x. When the periodic geodesies form a set of measure 
zero in S*M, the remainders are of order o(A'^~^), and one could average over the shorter 
interval [A, A + 1]. In the negatively curved case, 

fl(A,.)=0(^). 

Combining the Weyl and local Weyl law, we find the surface average of is a limit of traces: 

fiL{S*M) Js*M 

(43) ^ 

= i^iV(A) 

Xj<\ 

Here, diiL is Liouville measure on S*M ( §1.11) . 

3. Invariant states defined by eigenfuntions and their quantum limits 

When speaking of 'states' in quantum mechanics, one might refer to a normalized wave 
function or alternatively to the matrix elements p-^iA) = {Aip,ip) of an observable in the 
state. The latter use of 'state' is a standard notion in C* algebras, and is central in quantum 
ergodicity and mixing. The states evidently have the properties: (i) p^{A*A) > 0; {ii)p^{I) = 
1; {iii)p^ is continuous in the norm topology. The classical analogue is a probability measure, 
viewed as a positive normalized linear functional on C{S*M). In particular, eigenfunctions 
define states on (more correctly, its closure in the norm topology) as in f|T2l) . which we 
repeat: 

(44) pkiA) = {Aipk,ipk). 

We could (and will) also consider 'transition amplitudes' Pj,k{A) = {Aipj,ipk). 

It is an immediate consequence of the fact that U^^pj = e^^^^ipj that the diagonal states pk 
are invariant under the automorphism at f fT^ : 

(45) pk{UtAU;)=pu{A). 

In general, we denote by £ the compact, convex set of states in the vector space of continuous 
linear functionals on the closure of \E'° in its norm topology. We denote by 

(46) £u = {pe£:poat = p} 

the compact convex set of invariant states. For simplicity of notation, we continue to denote 
the closure by \Ef°, and refer to [Z6j for a detailed exposition. 

From a mathematical point of view, such the states Pk{A) are important because they 
provide the simplest means of studying eigenfunctions: one tests eigenfunctions against 
observables by studying the values pk{A). Many standard inequalities in PDE (e.g. Carleman 
estimates) have the form of testing eigenfunctions against well chosen observables. Of course, 
the high eigenvalue asymptotics is not simple. One would like to know the behavior as 
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Xj oo (or ^ of the diagonal matrix elements (Aipj, (pj) and the transition amplitudes 
{Aipi,ipj) between states for A G \E'°(M). One of the principal problems of quantum chaos 
is the following: 

Problem 3.1. Determine the set Q of 'quantum limits', i.e. weak* limit points of the 
sequence of invariant eigenfunction stats {pk} or equivalently of the Wigner distributions 
{Wk}. 

As will be illustrated below in simple examples, weak limits reflect the concentration and 
oscillation properties of eigenfunctions. 
Off-diagonal matrix elements 

(47) p,kiA) = {Aif^^ifj) 

are also important as transition amplitudes between states. They no longer define states 
since Pjk{I) = 0, are no longer positive, and are no longer invariant. Indeed, pjk{UtAU^) = 
Q^K^j-^k) p.^(^j\^^ SO they are eigenvectors of the automorphism at of ( |T^ . A sequence of such 
matrix elements cannot have a weak limit unless the spectral gap A-,- — \k tends to a limit 
r G M. Problem 13.11 has the following extension to off-diagonal elements: 

Problem 3.2. Determine the set Qr of 'quantum limits', i.e. weak* limit points of the 
sequence {Wkj} of distributions on the classical phase space S*M, defined by 

adWkj ■■= {Op{a)ipk, ipj) 

where Xj — Xk = t + o(l) and where a G C°°{S*M), or equivalently of the functional pjk. 

3.1. Simplest properties of quantum limits. The first is that {Kipki^k) for any 

compact operator K. Indeed, any orthonormal basis such as {^Pk} tends to zero weakly in L^. 
Hence {Kipk} tends to zero in for any compact operator and in particular the diagonal 
matrix elements tend to zero. It follows that for any A G \E'°(M), any limit of a sequence of 
{A^pk, fk) is equally a limit of {{A + K)(pk, fk)- 

Any two choices of Op, i.e. of quantizations of homogeneous symbols (of order zero) 
as pseudo-differential operators, are the same to leading order. Hence their difference is 
compact. Since a negative order pseudo-differential operator is compact, Q is independent 
of the definition of Op. 

These properties do not use the fact that (fj are eigenfunctions. The next property of the 
Pk are consequences of the fact pk E S^. 

Proposition 3.3. Q C Aii, where A4j is the convex set of invariant probability measures 
for the geodesic flow. They are also time-reversal invariant. 

To see this, we first observe that any weak * limit of of {pk} is a positive linear functional 
on symbols, which we identify with homogeneous functions of order zero on T*M or with 
smooth functions on S*M. Indeed, any limit of {Aipki^k)-, is bounded by inf^ ||A + 
(the infimum taken over compact operators), and for any A G \I'°, ||o"yi||Loo = inf;^ HA + i^'H. 
Hence any weak limit is bounded by a constant times ||crA||L°° and is therefore continuous 
on C{S*M). It is a positive functional since each pj is and hence any limit is a probability 
measure. The invariance under follows from pk G E^. by Egorov's theorem: any limit of 
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Pfc(A) is a limit of pk{Op{crA ° g^)) and hence the hmit measure is (?* invariant. Furthermore, 
the hmit measures are time-reversal invariant, i.e. invariant under (x, ^) {x, — ^) since the 
eigenfunctions are real-valued. 

Problem 13. II is thus to identify which invariant measures in A4j show up as weak limits of 
the functionals pk or equivalently of the Wigner distributions dWk- The problem is that Aij 
can be a very large set. Examples of invariant probability measures for the geodesic flow 
include: 

(1) Normalized Liouville measure dpL. In fact, the functional u of fl43p is also a state 
on \E^° for the reason explained above. A subsequence {v^j^.} of eigenfunctions is 
considered diffuse if pj^, — > u. 

(2) A periodic orbit measure p^ defined by p-yi^A) = J^aAds where is the length 
of 7. A sequence of eigenfunctions for which p^^ —>■ p^ obviously concentrates (or 
strongly 'scars') on the closed geodesic. 

(3) A finite convex combination YljLi Cjdp^^ of periodic orbit measures. 

(4) A mixed measure such as ^dpL + | J2jLi Cjdp-y. 

(5) A delta-function along an invariant Lagrangian manifold A C S*M. The associated 
eigenfunctions are viewed as localizing along A. 

(6) There are many additional kinds of singular measures in the Anosov case. 

Thus, the constraint in the Proposition is far from enough to determine Q. 

To pin down Q, it is necessary to find more constraints on the states pk and their weak* 
limits. We must use the quantum mechanics of p^ to pin down the possible classical limits. 

What possible additional constraints are there on the pfc? To date, only two are known. 

• Further symmetries of the ipj. In general there are none. But in special cases they 
exist. The most significant case is that of Hecke eigenfunctions, which carry an 
infinite number of further symmetries. See ^and |LINDl [Soundlj . 

• Entropies of limit measures. Lower bounds were obtained by N. Anantharaman and 
S. Nonnemacher ^ lANj in the case of {M,g) with Anosov geodesic flow (see also 
[XNKllRwj .) 

These constraints are of a very different nature. The Hecke symmetries are special to 
arithmetic hyperbolic quotients, where one supplement the geodesic flow-wave group con- 
nection with the Hecke symmetries and connections to number theory. The entropy bounds 
are very general and only use the geodesic flow-wave group connection. 

3.2. Ergodic sequences of eigenfunctions. A subsequence {fjf,} of eigenfunctions is 
called ergodic if the only weak * limit of the sequence of p^^ is dp^ or equivalently the 
Liouville state u. If dW^^ — > uj then in particular, we have 

1 fir M2,,r 1 Vol{E) 



Vol{M) Je ' Vol{M) 

for any measurable set E whose boundary has measure zero. In the interpretation of 
\ipk^{x)\'^dVol as the probability density of finding a particle of energy \\ at x, this says 
that the sequence of probabilities tends to uniform measure. 



20 



STEVE ZELDITCH 



However, Wk^ 00 is much stronger since it says that the eigenfunctions become diffuse 
on the energy surface S*M and not just on the configuration space M. One can quantize 
characteristic functions l^e of open sets in S*M whose boundaries have measure zero. Then 

{Op{lE)^Pj, ^Pj) = the amphtude that the particle in energy state hes in E. 

For an ergodic sequence of eigenfunctions, 

so that the particle becomes diffuse, i.e. uniformly distributed on S*M. This is the quantum 
analogue of the property of uniform distribution of typical geodesies of ergodic geodesic flows 
(Birkhoff's ergodic theorem). 

3.3. QUE. The Laplacian A or {M,g) is said to be QUE (quantum uniquely ergodic) if 
Q = {/iL}; i-e. the only quantum limit measure for any orthonormal basis of eigenfunctions 
is Liouville measure. 

Conjecture 3.4. (Rudnick-Sarnak, |RS] ) Let (M. g) be a negatively curved manifold. Then 
A IS QUE. 

In |RSj , the case of arithmetic manifolds is investigated and the role of the Hecke operators 
is clarified and exploited. In particular the arithmetic QUE Conjecture refers to the limits of 
Hecke eigen-states, i.e. eigenfunctions of the arithmetic symmetries called Hecke operators. 
As reviewed in ^ in this case the conjecture in all its forms is more or less completely 
solved. E. Lindenstrauss [LIND] (together with the recent final step in the noncompact 
case by Soundararajan |Soundlj ) proves this arithmetic QUE for arithmetic surfaces. If 
the multiplicity of the eigenvalues are uniformly bounded then one can deduce full QUE 
from arithmetic QUE, i.e. QUE for any orthonormal basis of eigenfunctions (see ^HD- Such 
uniform bounds on multiplicities are however far out of the range of current technology . 
The analogue of arithmetic QUE for Hecek holomorphic forms on noncompact arithmetic 
surfaces has recently been settled by Holowinsky and Soundararajan [Holj IHS] . In this case, 
the Hecke condition cannot be dropped due to the high multiplicity of such forms. Their 
methods are entirely arithmetical and we wont discuss them further here. 

Although we are not discussing quantum cat maps in detail, it should be emphasized that 
quantizations of hyperbolic (Anosov) symplectic maps of the torus) are not QUE. For a 
sparse sequence of Planck constants hk, there exist eigenfunctions of the quantum cat map 
which partly scar on a hyperbolic fixed point (see Faure-Nonnenmacher-de Bievre |FNBj ). 
The multiplicities of the corresponding eigenvalues are of order h'^^/\ iogh\. It is unknown if 
anything analogous can occur in the Riemannian setting, but as yet there is nothing to rule 
it out. 

As pointed out in |Z10j . QUE would follow if one could prove that quantum limits were 
invariant under a uniquely ergodic flow such as the horocycle flow of a compact hyperbolic 
quotient. A problem in trying to use this approach is that the horocycle flow is not Hamil- 
tonian with respect to the standard symplectic form on T*Xr, i.e. it cannot be quantized. 
It is however Hamiltonian with respect to a modified symplectic structure. The modified 
symplectic structure corresponds to letting the weight of automorphic forms vary with the 
A eigenvalue. As a result, one does get large sequences of automorphic forms which are 
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Liouville distributed. But QUE in the standard sense refers to sequences with fixed weight. 
It turns out that the weight has to grow so quickly for the QUE sequence that one cannot 
seem to relate the QUE result to quantum limits of forms with fixed weight [ZlOj . 

3.4. Simplest example: S^. The only computable example of ergodic eigenf unctions is 
the sequence of normalized eigenf unctions 2 sin kTix, 2 cos kipx on S^. Note that they are real 
valued; the exponentials e*'^^ are not quantum ergodic. 

The energy surface of T*S^ = S*^ x is the pair of circles |^| = ±1. The geodesic fiow has 
two invariant sets of positive measure (the two components), but the fiow is time- reversal 
invariant under (x,^) — > (x — ^) and the quotient fiow is ergodic. The real eigenfunctions are 
invariant under complex conjugation and therefore the quotient quantum system is ergodic. 

Quantum ergodicity in this case amounts to 

— / V{x){smkxYdx ^ — / Vdx. 
Jo 27r Jo 

This is obvious by writing sin kx in terms of exponentials and using the Riemann-Lebesgue 
Lemma. 

This simple example illustrates an important aspect of quantum limits: They are weak lim- 
its which owe to the fast oscillation of the eigenfunction squares. In the limit the oscillating 
functions tend to their mean values in the weak sense. 

In particular, it illustrates why we only consider squares and not other powers of eigen- 
functions: weak limits are not preserved under non-linear functionals such as powers. The 
study of LP norms of chaotic eigenfunctions is very difficult (there are Iwaniec-Sarnak for 
Hecke eigenfunctions). 

3.5. Example of quantum limits: the fiat torus. The only examples where one can 
compute quantum limits directly are the completely integrable ones such as the standard 
sphere, torus or symmetric spaces. These examples of course lie at the opposite extreme from 
chaotic or ergodic dynamics. We use the simplest one, the fiat torus ]R"/Z", to illustrate the 
definition of weak* limits. 

An orthonormal basis of eigenfunctions is furnished by the standard exponentials e^'^*^'^'^^ 
with k E Z". Obviously, |e^'^*^'^'^^ p = 1, so the eigenfunctions are diffuse in configuration 
space. But they are far from diffuse in phase space. For any pseudodifferential operator, 
j^^2m(k,x) _ ^^^^ j^-^^2m{k,x) -yy-j-^grg (i(^x^ k) is the complete symbol. Thus, 

^^g2.i(M>^e2.i(M>^ = /■ a{x,k)dx^ [ aA{x,^)dx. 

Thus, the Wigner distribution is A subsequence e^'^*^^^'^^ of eigenfunctions has a weak 
limit if and only if tends to a limit vector in the unit sphere in M". In this case, the 
associated weak* limit is f^„^^„aA{x,^o)dx, i.e. the delta-function on the invariant torus 
T^Q G S*M defined by the constant momentum condition ^ = ^q. The eigenfunctions are 
said to localize under this invariant torus for g^. 

The invariant torus is a Lagrangian submanifold of T*M"/Z", i.e. a submanifold of dimen- 
sion n on which the standard symplectic form dx A d^ restricts to zero. The exponentials 
are special cases of WKB or Lagrangian states a/i.e^'^, where is a semi-classical symbol 
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an, ~ ajh^ . The associated Lagrangian submanifold is the graph (x, dS{x)) of S. Thus, {k, x) 
generates the Lagrangian submanifold ^ = jfj- 

In general, one says that a sequence {^j^} of eigenfunctions concentrates microlocally on 
a Lagrangian submanifold K (Z S*M ii 

(48) {A^j^,^j^) / aAdv, 

J A 

for some probability measure u on A. Necessarily, A is invariant under g^. 

3.6. Hyperbolic case. Using the Helgason Poisson integral formula, the Wigner distribu- 
tions can be expressed in terms of the dTir^ and Cir^ as follows. As in [Z3j . we define the 
hyperbolic calculus of pseudo-differential operators Op{a) on D by 

Op(a)e(5+-)<^.^) = a(z,6,r)e(5+"')(-,f'>. 

We assume that the complete symbol a is a polyhomogeneous function of r in the classical 
sense that 

oo 

a{z,b,r) ~ J]aj(z,6)r-^+™ 
i=o 

for some m (called its order). By asymptotics is meant that 

R 

a{z, b,r)-J2 «i(^. b)r-^+"' G 

j=0 

where a G S''^ if snp{l + ry~''\D"D^Dia{z, 6, r)| < +oo for all compact set and for all a, P,j. 

The non-Euclidean Fourier inversion formula then extends the definition of Op{a) to 
C-(D): 

Op{a)u{z) = / / a{z,b,r)e^ 2+-'^){^^^) J^u{r, 6)rtanh(27rr)dr|d6|. 

J B Jr 

A key property of Op is that Op{a) commutes with the action of an element 7 G G 
(T^m(z) = ^(7-2)) if and only if a{jz,'yb,r) = a{z,b,r). F-equivariant pseudodifferential 
operators then define operators on the quotient Xp. 

By the Helgason-Poisson formula one has a relative explicit formula for the Wigner dis- 
tributions Wir. G V'{S*Xr)) defined by 

(49) (a, W„,) = I a{g)W„.M9) ■= {Op{a)<f„„ f^r,) m^Lr), « e C°°(5*Xr). 

Js*Xr 

Equivalently we have 

(50) Wir^ = (fiir^eiry 

Thus, the Wigner distributions are far more diffuse in phase space than in the case of the flat 
torus. But this does not rule out that their weak* limits could have a singular concentration. 
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3.6.1. Patterson- Sullivan distributions. Egorov's theorem implies that Wigner distributions 
tend to invariant measures for the geodesic flow. The question arises whether there exist g^- 
invariant distributions constructed from eigenfunctions which are asymptotic to the Wigner 
distributions. In [AZj such distributions were constructed on hyperbolic surfaces and termed 
"Patterson- Sullivan distributions" by analogy with their construction of boundary measures 
associated to ground states on infinite volume hyperbolic manifolds. 

Definition: The Patterson-Sullivan distribution associated to a real eigenfunction ipij.. is 
the distribution on B x B \A defined by 

^ , ^ Tir.[dh)Tir\dh') 

pSir,{db , dh) := ^^_y^r+2ir, 



If ifir^ is F-automorphic, then psirj{db' ,db) is P-invariant and time reversal invariant. 

To obtain a 5f*-invariant distribution on S'*Xr, we tensor psirj with dt. We then normalize 
by dividing by the integral against 1. The result is an invariant distribution PSir^ for 
constructed as a quadratic expression in the eigenfunctions. In |AZ] Theorem 1.1, it is proved 
(theorem that 

f aPS„, = I aW,r,+0{r-^). 

J 5'*Xp J iS'^'Xp 

Hence the quantum limits problem is equally one of determining the weak* limits of the 
Patterson-Sullivan distributions. It is shown in |AZ] that they are residues of dynamical 
L-functions and hence have a purely classical definition. 

In fact, there is an explicit intertwining operator L^.^ mapping PSi^^ ^irj have 

(51) (a, W„^) = (a, pX,) + rJ^L^ia), PS,,.,) + 0{rf). 

4. Quantum ergodicity and mixing of eigenfunctions 

In this section, we review a basic result on quantum ergodicity. We assume that the 
geodesic flow of (M, g) is ergodic. Ergodicity of g^ means that Liouville measure dfj,L is an 
ergodic measure for (7* on S*M, i.e. an extreme point of Aij. That is, there any (^'-invariant 
set has Liouville measure zero or one. Ergodicity is a spectral property of the operator V^* of 
(IT^ on L'^{S*M, dfii), namely that V* has 1 as an eigenvalue of multiplicity one. That is, the 
only invariant functions (with respect to Liouville measure) are the constant functions. 

In this case, there is a general result which originated in the work of A. I. Schnirelman and 
was developed into the following theorem by S. Zelditch, Y. Colin de Verdiere on manifolds 
without boundary and by P. Gerard-E. Leichtnam and S. Zelditch-M. Zworski on manifolds 
with boundary. 

Theorem 4.1. Let {M,g) be a compact Riemannian manifold (possibly with boundary), and 
let {Xj^Lpj} be the spectral data of its Laplacian A. Then the geodesic flow G* is ergodic on 
{S*M, dni) if and only if, for every A G \1/°(M), we have: 

• (l) limA^oojvPo EA,<Al(^<^i'<^i) = 0- 

• (^^) (Ve)(35)limsupA_^^^,^fc:A,,Afe<A \{A^j,^k)? < e 
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The diagonal result may be interpreted as a variance result for the local Weyl law. Since all 
the terms are positive, the asymptotic is equivalent to the existence of a s subsequence {^j^} 
of eigenfunctions whose indices jk have counting density one for which {Atpj^^^pj^) — >• uj{A) 
for any A G \E'°(M). As above, such a sequence of eigenfunctions is called ergodic. One can 
sharpen the results by averaging over eigenvalues in the shorter interval [A, A + 1] rather than 
in [0,A]. 

The off-diagonal statement was proved in |Z9] and the fact that its proof can be reversed 
to prove the converse direction was observed by Sunada in [Suj . A generalization to finite 
area hyperbolic surfaces is in |Z8j . 

The first statement (i) is essentially a convexity result. It remains true if one replaces the 
square by any convex function F on the spectrum of A, 

(52) ^(^^V^'^' V'^) - ^(^)) ^ 0- 

Xj<E 



4.1. Quantum ergodicity in terms of operator time and space averages. The diag- 
onal variance asymptotics may be interpreted as a relation between operator time and space 
averages. 

Definition Let A E'^^ he an observable and define its time average to be: 

{A) := hm {A)t. 

I — ^oo 



where 

U'AU-'dt 



and its space average to be scalar operator 

uj{A) ■ I 

Then Theorem 14.11 (1) is (almost) equivalent to, 

(53) {A) = uj{A)I + K, where lim uJx{K*K) 0, 

A— >oo 

where uj\{A) = jjj^TrE{\)A. Thus, the time average equals the space average plus a term K 

which is semi-classically small in the sense that its Hilbert-Schmidt norm square ||£'A-ft'|||^5 
in the span of the eigenfunctions of eigenvalue < A is o(iV(A)). 

This is not exactly equivalent to Theorem 14.11 (1) since it is independent of the choice of 
orthonormal basis, while the previous result depends on the choice of basis. However, when 
all eigenvalues have multiplicity one, then the two are equivalent. To see the equivalence, note 
that {A) commutes with a/A and hence is diagonal in the basis {^fj} of joint eigenfunctions 
of (A) and of Ut- Hence K is the diagonal matrix with entries {Aipk,(pk) — uj{A). The 
condition is therefore equivalent to 

^ ^ \i<E 
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4.2. Heuristic proof of Theorem 14.11 (i) . There is a simple picture of eigenfunction states 
which makes Theorem 14.11 seem obvious. Justifying the picture is more difficult than the 
formal proof below but the reader may find it illuminating and convincing. 

First, one should re- formulate the ergodicity of (7* as a property of the Liouville measure 
d^L- ergodicity is equivalent to the statement dfi^ is an extreme point of the compact convex 
set Aij. Moreover, it implies that the Liouville state u on \1/°(M) is an extreme point of the 
compact convex set S^. of invariant states for at of (fT4l) : see ^uj for background. But the 
local Weyl law says that uj is also the limit of the convex combination tv^e) Ylx <EPj- 
extreme point cannot be written as a convex combination of other states unless all the states 
in the combination are equal to it. In our case, uj is only a limit of convex combinations so 
it need not (and does not) equal each term. However, almost all terms in the sequence must 
tend to UJ, and that is equivalent to (1). 

One could make this argument rigorous by considering whether Liouville measure is an 
exposed point of Sj and Aij. Namely, is there a linear functional A which is equal to zero 
at UJ and is < everywhere else on £f7 If so, the fact that j^^J2\j<E-^(Pj) ~^ ^ implies 
that A(pj) ^ for a subsequence of density one. For one gets an obvious contradiction if 
^{Pjk) < — e < for some e > and a subsequence of positive density. But then pj^, uj 
since uj is the unique state with A(p) =0. 

In |J] it is proved that Liouville measure (or any ergodic measure) is exposed ui M.j. It is 
stated in the following form: For any ergodic invariant probability measure /i, there exists 
a continuous function f on S*M so that fi is the unique /-maximizing measure in the sense 
that 

J fdjji = sup I J fdm : m G Aij 

To complete the proof, one would need to show that the extreme point uj is exposed in Sj 
for the C* algebra defined by the norm-closure of \1/°(M). 

4.3. Sketch of Proof of Theorem 14.11 (i). We now sketch the proof of fl52|) . By time 
averaging, we have 

(54) ^F((A^,,^,)-cu(A))= ^ F(((A)T-cu(A)y;fe,^,)). 

Xj<E \j<E 

We then apply the Peierls-Bogoliubov inequality 

n 

Y,F{{B^„^,))<TtF{B) 
with B = Ue[{A)t - uj{A)]Ue to get: 

(55) Yl - '^{A)v>k, fk)) < Tr F{Ue[{A)t - uj{A)]Ue). 

\,<E 

Here, H^; is the spectral projection for H corresponding to the interval [0,-E']. By the 
Berezin inequality (if -F(O) = 0): 

j^i:vF{YiE[{A)T-uj{A)]YlE) < j^TvTiEF{[{A)T-uj{A)])YiE 



= uje{v{{A)t-oj{A))). 
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As long as F is smooth, F{{A)t — ^{A)) is a pseudodifferential operator of order zero with 
principal symbol F{{aA)T — ^^{A)). By the assumption that ue ^ uj we get 

We) S ~ ""^^^^ - / ^^^^^^^ ~ 

^ ^ \,<E •J\H=1} 



As T — s> oo the right side approaches <y9(0) by the dominated convergence theorem and by 
Birkhoff 's ergodic theorem. Since the left hand side is independent of T, this implies that 



J"^ WW) S ^^^^^^^ '^^^ ' ^(^)) = 



for any smooth convex F on Spec(y4) with -F(O) = 0. □ 

This proof can only be used directly for scalar Laplacians on manifolds without boundary, 
but it still works as a template in more involved situations. For instance, on manifolds 
with boundary, conjugation by the wave group is not a true automorphism of the observable 
algebra. In quantum ergodic restriction theorems (see the appropriate conjugation is an 
endomorphism but not an automorphism. Or when A has continuous spectrum (as in finite 
area hyperbolic surfaces), one must adapt the proof to states which are not L^-normalized 



4.4. QUE in terms of time and space averages. The quantum unique ergodicity prob- 
lem (the term is due to Rudnick-Sarnak |RSj ) is the following: 

Problem 4.2. Suppose the geodesic flow g*' of{M,g) is ergodic on S*M. Is the operator K 
in 

{A) = uj{A) + K 

a compact operator? Equivalently is Q = {dfii}? In this case, ^/A is said to be QUE 
(quantum uniquely ergodic) 

Compactness of K implies that {Kifk, fk) 0, hence {A(pk, fk) ~^ ^{A) along the entire 
sequence. 

Rudnick-Sarnak conjectured that A of negatively curved manifolds are QUE, i.e. that for 
any orthonormal basis of eigenfunctions, the Liouville measure is the only quantum limit 

m- 

4.5. Converse QE. So far we have not mentioned Theorem 14.11 (2). An interesting open 
problem is the extent to which (2) is actually necessary for the equivalence to classical 
ergodicity. 



Problem 4.3. Suppose that y/A is quantum ergodic in the sense that (1) holds in Theorem 
4.1\ What are the properties of the geodesic flow g^. Is it ergodic (in the generic case)? 



In the larger class of Schrodinger operators, there is a simple example of a Hamiltonian 
system which is quantum ergodic but not classically ergodic: namely, a Schrodinger operator 
with a symmetric double well potential W. That is, is a shaped potential with two 
wells and a Z2 symmetry exchanging the wells. The low energy levels consist of two connected 
components interchanged by the symmetry, and hence the classical Hamiltonian flow is not 
ergodic. However, all eigenfunctions of the Schrodinger operator — + W are either even 
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or odd and thus have the same mass in both wells. It is easy to see that the quantum 
Hamiltonian is quantum ergodic. 

Recently, B. Gutkin [GutJ has given a two dimensional example of a domain with boundary 
which is quantum ergodic but not classically ergodic and which is a two dimensional analogue 
of a double well potential. The domain is a so-called hippodrome (race-track) stadium. 
Similarly to the double well potential, there are two invariant sets interchanged by a Z2 
symmetry. They correspond to the two orientations with which the race could occur. Hence 
the classical billiard flow on the domain is not ergodic. After dividing by the J? symmetry 
the hippodrome has ergodic billiards, hence by Theorem l4.H the quotient domain is quantum 
ergodic. But the The eigenf unctions are again either even or odd. Hence the hippodrome is 
quantum ergodic but not classically ergodic. 

Little is known about converse quantum ergodicity in the abscence of symmetry. It is 
known that if there exists an open set in S*M filled by periodic orbits, then the Laplacian 
cannot be quantum ergodic (see [MOZj for recent results and references). But it is not even 
known at this time whether KAM systems, which have Cantor-like invariant sets of positive 
measure, are not quantum ergodic. It is known that there exist a positive proportion of 
approximate eigenfunctions (quasi-modes) which localize on the invariant tori, but it has 
not been proved that a positive proportion of actual eigenfunctions have this localization 
property. 

4.6. Quantum weak mixing. There are parallel results on quantizations of weak-mixing 
geodesic flows. We recall that the geodesic flow of (M, g) is weak mixing if the operator 
has purely continuous spectrum on the orthogonal complement of the constant functions in 
L\S*M, d^L)- 

Theorem 4.4. |Z8j The geodesic flow of (M, g) is weak mixing if and only if the conditions 
(l)-(2) of Theorem \4.1\ hold and additionally, for any A G ^°(Af), 

(Ve)(35)limsup-i- \{A^,,^,)\' < e (Vr G M) 

The restriction j ^ k is of course redundant unless r = 0, in which case the statement 
coincides with quantum ergodicity. This result follows from the general asymptotic formula, 
valid for any compact Riemannian manifold {M,g), that 



(56) 



sinT(Ai 



r(A,-A,-r) 



^ f_^e^'^V,{aA)\\l-\^\MAf 



I 2T 

In the case of weak-mixing geodesic flows, the right hand side — > as T ^ 00. As with 
diagonal sums, the sharper result is true where one averages over the short intervals [A, A + 1]. 

Theorem 14.41 is based on expressing the spectral measures of the geodesic flow in terms of 
matrix elements. The main limit formula is: 
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where rf/io-^ is the spectral measure for the geodesic flow corresponding to the principal 
symbol of A, a a ^ C°^(S'*M, d[ii). Recall that the spectral measure of V* corresponding to 
f E L'^ is the measure dfif defined by 

Jr 

4.7. Evolution of Lagrangian states. In this section, we briefly review results on evo- 
lution of Lagrangian states and coherent states. We follow in particular the article of R. 
Schubert |Schu3j . 

A simple Lagrangian or WKB state has the form = a{h, x)eTi^^^^ where a{h,x) is a 

semi-classical symbol a ~ JZ'jLo fi^^ciji^x). The phase S generates the Lagrangian submanifold 
{x,dS{x)) ^T*M. 

It is proved in Theorem 1 of |Schu3j that if is Anosov and if A is transversal to the 
stable foliation (except on a set of codimension one), then there exists C, r > so that 
for every smooth density on A and every smooth function a G C°°{S*M), the Lagrangian 
state tp with symbol satisfies. 



(5J 



S*M J A 



In order that the right side tends to zero as ^ — > 0, t — > cxd it is necessary and sufficient that 

t<^-^\\ogh\. 



5. Concentration of eigenfunctions around hyperbolic closed geodesics 

As mentioned above, the quantum ergodicity Theorem 14.11 leaves open the possible ex- 
istence of a sparse (zero density) subsequence of eigenfunctions which 'weakly scar' on a 
hyperbolic orbit 7 in the sense that its quantum limit z/q contains a non-trivial multiple 
of the periodic orbit measure c/i^ as a non-zero ergodic component. The Anantharaman 
entropy bound shows that when (M, g) is Anosov, there cannot exist such a sequence of 
eigenfuntions (or even quasi-modes) which tend to yU^ itself, but a quantum limit could have 
the form ci/i^ + C2/iL for certain ci, C2 satisfying ci + C2 = 1. The question we address in this 
section is the possible mass profile of such a scarring eigenfunction in a neighborhood of 7. 
More precisely, how much mass does (fx have in a shrinking h}^"^""^ tube around a hyperbolic 
closed geodesic? This question will surface again in §1L21[ 

We first note that the existence of a quantum limit of the form ci/i^ + 02^1 for {M,g) 
with Anosov geodesi flow is not so implausible. Analogous eigenfunction sequences do exist 
for the so-called quantum cat map |FNBl [FN]. And exceptional sequences 'scarring' on a 
certain 1-parameter family of periodic orbits exists for the Bunimovich stadium [Has] . At 
this time, there are no known examples of sequences of eigenfunctions or quasi-modes for 
(M, g) with ergodic geodesic flow that 'weakly scar' along a hyperbolic closed orbit 7 and 
no results prohibiting them. 

For simplicity, assume that (M, g) is a Riemannian manifold of dimension 2 with a closed 
geodesic 7 of length L. We assume 7 is an embedded (non self-intersecting) curve. If (pj is 
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a Laplace eigenf unction, we define the mass profile of ^pj near 7 to be the function 

J d(x,'y)=r 

where dS A dr = dV. 

Before considering possible mass profiles of eigenfunctions scarring on hyperbolic closed 
geodesies, let us recall the opposite and much better known case of scarring of Gaussian 
beams along elliptic closed geodesies on surfaces [Ral IRa2l IB.B] . When 7 is an elliptic 
closed geodesic, then there always exists a sequence of quasi-modes (Gaussian beams) which 
concentrates on 7. As the name suggests, Gaussian beams tpx oscillate like e*'*'* along 7 
and resemble Gaussians \/Ae~^^"^*^*-'^'^^ (for some positive symmetric matrix A{s)) in the 
transverse direction with height A^/^ and concentrated in a tube around 7. Thus, the 

mass profile is a Gaussian probability measure with mean zero and variance A~^/^. The local 
model for such quasi-modes is that of a harmonic oscillator in the fibers of A^^. One can 
construct the quasi-mode so that it is of infinite order. Of course, stable elliptic orbits do 
not exist when {M,g) has ergodic geodesic flow (by the KAM theorem). 

Now consider hyperbolic closed orbits. It was pointed out by Duistermaat [Duij (sec- 
tion 1.5) that one cannot construct analogous quasi-modes associated to hyperbohc closed 
geodesies as Lagrangian states. The stable/unstable manifolds A± of 7, containing the 
geodesies which spiral in towards 7, are invariant Lagrangian submanifolds, but the only 
invariant half-density on the Lagrangians is the 'delta'-density on the closed geodesic. 

Further, there are apriori limitations on the degree to which eigenfunctions sequences can 
concentrate around hyperbolic closed geodesic on any {M,g). 

Theorem 5.1. \BZ\ IChrt IGVPj Let {M,g) be a compact Riemannian manifold, and let 7 
he a hyperbolic closed geodesic. Let U be any tubular neighborhood of 'j in M. Then for any 
eigenfunction ipx, there exists a constant C depending only on U such that 

[ \vx\'dV,>-^\\y^,\\l,. 
Jm\u log A 

More generally, let A G \1/°(M) be a pseudo-differential orbit whose symbol equals one in 
a neighborhood of 7 in S*M and equals zero outside another neighborhood. Then for any 

eigenfunction (^x - A)lpx\\l2 > 7§il Iv^aI |l2- 

This allows scarring sequences along a hyperbolic orbit to occur, it just limits the rate 
at which the mass concentrates near 7. It implies that the mass profile of a sequence of 
eigenfunctions concentrating on a hyperbolic close geodesic has "long tails", i.e. there is a 
fairly large amount of mass far away from the geodesic, although sequences with the quantum 
limit fi.y must tend to zero outside of any tube around the closed geodesic. Note that this 
result makes no dynamical hypotheses. It applies equally to (M, g) with integrable geodesic 
flow and to those with Anosov geodesic flow. To the author's knowledge, there do not exist 
more precise results in the Anosov case. 

An obvious question at this point is whether there are any examples of (M, g) , with any 
type of geodesic flow, possessing a sequence of eigenfunctions scarring on a closed hyperbolic 
orbit. The answer to this question is 'yes'. It is simple to see that such eigenfunctions exist 
in the opposite extreme of completely integrable systems, for instances surfaces of revolution 
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like peanuts with hyperbolic waists. Examples include joint eigenf unctions of the square root 
of the Laplacian and rotation on surfaces of revolution with a hyperbolic waist. A truncated 
hyperbolic cylinder is another example studied in |CVPj . In this case, the joint spectrum 
fills out the image of the moment map {pg, |^|) : T*M — >• M^, where pe{x.,i) = {i, ^) is the 
angular momentum. At critical distances to the axis of rotation, the lattitude circles are 
closed geodesies and the level set of the moment map becomes singular. If the surface is 
shaped like a peanut, the waist is a hyperbolic closed geodesic. Joint eigenfunctions whose 
joint eigenvalues are asymptotic to the singular levels of the moment always exist. The 
modes concentrate on the level sets of the moment map, and in fact they concentrate on the 
hyperbolic closed geodesic. 

5.1. Mass concentration of special eigenfunctions on hyperbolic orbits in the 
quantum integrable case. The mass profile of scarring eigenfunctions near a hyperbolic 
in the completely integrable case is studied in jCVP] on tubes of fixed radius and in |NVIITZ2] 
on tubes of shrinking radius. Let 7 C S*M be a closed hyperbolic geodesic of an (M, g) with 
completely integrable geodesic flow and for which is quantum integrable (i.e. commutes 
with a maximal set of pseudo-differential operators; see [TZ2j for background). We then 
consider joint eigenfunctions Ag and of these operators. It is known (see e.g. |TZ2j . Lemma 
6) that there exists a special sequence of eigenfunctions concentrating on the momentum 
level set of 7. We will call them (in these notes) the '-^-sequence. 

Assume for simplicity that the moment level set of 7 just consists of the orbit together with 
its stable/unstable manifolds. Then it is proved in |TZ2] that the mass of (p^ in the shrinking 
tube of radius around 7 with 5 < | is ~ (1 — 2(5) (see also [NVj for a closely related result in 
two dimensions). Thus, the mass profile of such scarring integrable eigenfunctions only differs 
by the numerical factor (1 — 25) from the mass profile of Gaussian beams. The difference 
is that the 'tails' in the hyperbolic case are longer. Also the peak is logarithmically smaller 
than in the elliptic somewhat weaker statement is proved in |TZ2j ). 

Let us state the result precisely and briefly sketch the argument. It makes an interesting 
comparison to the situation discussed later on of possible scarring in the Anosov case. 

We denote by vr : S*M M the standard projection and let vr(7) be the image of 7 in 
M. We denote by Te(7r(A)) the tube of radius e around 7r(A). For < 5 < 1/2, we introduce 
a cutoff Xi(a;; h) e C^{M) with < < 1, satisfying 

• (i) supp Xi C T;j^(7r(7)) 

• (ii) X? = 1 on T3/4^^(7r(7)). 

Theorem 5.2. Let be a hyperbolic closed orbit in {M,g) with quantum integrable Ag, and 
let {(f^} be an L? normalized '^-sequence of joint eigenfunctions Then for any < 5 < 1/2, 
\im.n,-^Q{Opn{x{)Vti^V^i) > (1 - 25). 

5.1.1. Outline of proof. For simplicity we assume dimM = 2. Let xH^i^'^ ^) ^ C^(T*M; [0, 1]) 
be a second cutoff supported in a radius tube, Q{h), around 7 with Q{h) C suppxi and 
such that Xi = 1 oil suppxi- Thus, xii^^O — xii^^O^ for any (x, ^) G T*M. By the 
Garding inequality, there exists a constant Ci > such that: 

(59) iOpnixi)^,, ^,)) > (Opnixi)^,, V,) - C,h'-^' . 

We now conjugate the right side to the model setting of 5*^ x M^, i.e. the normal bundle 
to 7. The conjugation is done by h Fourier integral operators and is known as conjugation 
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to quantum Birkhoff normal form. In the model space, the conjugate of is a function 
oi Ds = ^ along and the dilation operator := h{Dyy + yDy) along M. By Egorov's 
theorem 



(60) 



C.h 



1-25 



where u^{y, s; K) is the model joint eigenfunction of Dg, J'^, and c{K) is a normalizing constant. 
This reduces our problem to estimating the explicit matrix elements (Op;i(x2°'^)w^, u^) of the 
special eigenfunctions in the model setting. The operator has a continuous spectrum with 
generalized eigenfunctions y-i/2+»^/ft ^j^e eigenfunctions on the 'singular' level 7 correspond 
to A ~ Eh. A caluclation shows that the mass in the model setting is given by 



(61) 



log h 



X 



-y^+'>^/\[x/h^i)dx 



Classical analysis shows that the right side tends to 1 — 25 as ^ — if A ~ Eh. 



5.2. Comparison to Anosov case. This large mass profile may be a special feature of 
integrable systems, reflecting the fact that the stable and unstable manifolds of the hyperbolic 
closed orbit coincide. A heuristic picture of the mass concentration in this case is as follows: 
Since the eigenfunction is a stationary state, its mass must be asymptotically invariant under 
the geodesic flow. Since the flow compresses things exponentially in the stable direction and 
expands things in the unstable direction, the mass can only concentrate on the fixed closed 
geodesic and on the unstable manifold W^. But returns to 7 as the stable manifold 
(like a figure 8). Hence, the only invariant measure is the one supported on the closed 
geodesic and the the mass can only concentrate there. Although it does not seem to have 
been proved in detail yet, it is very plausible that the mass concentration in the integrable 
case provides an upper bound for any {M,g), i.e. it has 'extremal' mass concentration. 

In the Anosov case, the stable and unstable manifolds of 7 are transverse, so the dynamical 
picture is completely different. First, there is no obvious mechanism as in the integrable case 
forcing mass of any sequence of eigenfunctions to concentrate on the Lagrangian manifold 
formed by 7 and in the Anosov case. If mass did concentrate around 7, it would still be 
forced to concentrate on 7 and on W^, but becomes dense in S*M. Hence some of the 
mass must spread out uniformly over S*M and is lost from a neighborhood of 7. This makes 
it plausible that one does not get mass concentration for eigenfunctions around hyperbolic 
closed orbits of Anosov systems. 

Yet, in the "cat map" analogue, there do exist scarring eigenfunctions {FNBt IFN] for a 
special sparse sequence of Planck constants. The multiplicities of the eigenvalues of the 
cat map for this sequence saturate the bound h^^ /\\ogh\, and therefore one can build up 
eigenfunctions with very unexpected properties. There is a surprising quantum mechanism 
forcing producing concentration of special modes at hyperbolic fixed points which was dis- 
covered by Faure-Nonnenmacher-de Bievre. The eigenfunction amplitude spreads out along 
h^^^'^ segments of W^. These segments are h^^'^ dense and so there is interference between 
the amplitudes on close pieces of W^. The interference is constructive along W''' and the 
mass builds up there and then as in the integrable case gets recycled back to the hyperbolic 
orbit. 
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It is not known whether this phenomenon occurs in the Riemannian setting. It is pre- 
sumably related to the existence of sparse subsequences of eigenvalues with the same large 
multiplicities. 

6. Boundary quantum ergodicity and quantum ergodic restriction 

In this section, we briefly discuss quantum ergodic restriction theorems. The general 
question is whether restrictions of quantum ergodic eigenf unctions to hypersurfaces (or mi- 
crolocally, to cross-sections of the geodesic flow) are quantum ergodic on the hypersurfaces. 
We only consider here the case where the hypersurface is the boundary of a domain with 
boundary, which was studied in \GL\ IHZt IBuj . More general quantum ergodic restriction 
theorems are given in |TZ3j . The boundary results play an important role in the recent 
scarring results of A. Hassell for eigenfunctions on the stadium. 

We thus consider the boundary values of interior eigenfunctions 

A Uj = A| Uj in Q, {uj,Uk)L^n) = Sjk, 
Buj\Y = 0, Y = dn 

of the Euclidean Laplacian A on a compact piecewise smooth domain C M" and with 
classically ergodic billiard map /3 : B*Y B*Y, where Y = dfl. Here is a zeroth order 
semiclassical pseudo differential operator on Y. The relevant notion of boundary values 
(i.e. Cauchy data) depends on the boundary condition. We only consider the boundary 
conditions 

u\y, Dirichlet 

(62) Bu 



d^ulY, Neumann 



Let Ab denote the positive Laplacian on Q with boundary conditions Bu = 0. Then A^ 
has discrete spectrum 0<Ai<A2<---^cxd, where we repeat each eigenvalue according 
to its multiplicity, and for each Xj we may choose an normalized eigenfunction uj. 

The algebra of observables in the boundary setting is the algebra ^l/°(^) of zeroth order 
semiclassical pseudodifferential operators on Y, depending on the parameter h G [0, ho]. We 
denote the symbol of A = Ah E "^liY) by a = a{y,ri,h). Thus a{y,ri) = a{y,ri,0) is a 
smooth function on T*Y. 

To each boundary condition B corresponds 

• A specific notion of boundary value u'j of the eigenfunctions Uj. We denote the 
L^-normalized boundary values by -u^ = -u^/IIm^H. 

• A specific measure dfiB on B*{Y). 

• A specific state ujb on the space "^^{Y) of semiclassical pseuodifferential operators of 
order zero defined by 

(63) = [ a{y,7])dfiB. 

vo1(d" ^)vol(i2) Jb*y 

Here is a table of the relevant boundary value notions. In the table, da is the natural 
symplectic volume measure on B*Y. We also define the function 7(g) on B*Y by 

(64) 7(g) = - \r]\^, q = {y,r]). 
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Boundary Values 


B 


Bu 




dfiB 


Dirichlet 


u 


Y 




-f{q)da 


Neumann 




Y 


U Y 





The limit states are determined by dictated by the local Weyl law for the boundary 
condition B. 

Lemma 6.1. Let he either the identity operator on Y or a zeroth order semiclassical 
operator on Y with kernel supported away from the singular set. Then for any of the ahove 
houndary conditions B, we have: 

lim ^rTTTT {AhjU^j,u^j) = u!b{A), B = Neumann, 



(65) 



A-oo N{X) 



A,<A 



lim — — y {AhV., u)) =ub{A), B= Dirichlet. 



A,<A 



The boundary quantum ergodic theorem is: 



Theorem 6.2. [HZl ED ^ Let n c 



he a hounded piecewise smooth manifold with 



on 



lim {Ah^u^j,u^j) = iOsiA), B= Dirichlet, 



ergodic hilliard map. Let {Uj} he the houndary values of the eigenfunctions {uj} of Ab 
L'^{Q) in the sense of the tahle ahove. Let Ah he a semiclassical operator of order zero on Y . 
Then there is a suhset S of the positive integers, of density one, such that 

lim {AhM^j,u^j) = idBiA), B= Neumann, 

(66) 

where hj = Xj^ and ujb is as in fl63l) . 

In the case A = I and for the Neumann boundary condition, we have 

while for the Dirichlet boundary condition, 



'jWLHY) 



lim 

j-*oo,j£S 



II L2(y) 



vol(fi) ' 
2vol(F) 



6.1. Sketch of proof. The fact that the quantum limit state ub and the corresponding 
measure dfiB are not the natural symplectic volume measure da on B*Y is due to the fact 
that the quantum dynamics is defined by an endomorphism rather than an automorphism of 
the observable algebra. In the Neumann case, the dynamics are generated by the operator 
Fh on Y with kernel 

Fh{y,y') = 2^Go{y,y',h-'), y^y'eY, where 



(67) 



dl^y 

Go{y,y',X) = '-X''-\2nX\z - z'\p'^-'y'}ie^^;i^_^{X\z - z'\) 
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is the free outgoing Green function on R'^. These are (almost) semi-classical Fourier integral 
operators whose phase functions, the boundary distance function db{y, y') = \y—y' {y, y' G Y) 
generates the billiard map. It is well known that this operator leaves the boundary values 
of Neumann eigenf unctions invariant: 

(68) Fh,u'; = u'v J = 1'2,... 
It follows that the states 

(69) p,{A) := {A,m'^,u'^) 

are invariant for Fhj. The family {Fh} defines a semiclassical Fourier integral operator 
associated to the billiard map f3 (for convex Q), plus some negligible terms. The quantum 
dynamics on "^^(Y) is thus generated by the conjugation 

(70) a,^{A,^) = F*^ A,^ F,. 

This is analogous to the interior dynamics generated by the wave group U*^ with the boundary 
conditions, but unlike Z7|j, F^ is not unitary or even normal. Indeed, the zeroth order part 
of F^Fh is a semiclassical pseudodifferential operator with a non-constant symbol. 

The Egorov type result for the operator Fh is as follows: Let (3 denote the billiard map 
on B*Y° and let Ah = Op{ah) be a zeroth order operator whose symbol a{y, r], 0). Let 7 be 
given by flM|) . Then 

FhAhFh = Ah + Sh, 

where Ah is a zeroth order pseudodifferential operator and ||S'ft||L2^L2 < Ch. The symbol of 
Ah is 



(71) 



l{qh{(3{q))-Hm), q e B*Y 
0, B*Y. 



This is a rigorous version of the statement that Fh quantizes the billiard ball map. This 
Egorov theorem is relevant to the Neumann boundary problem. In the Dirichlet case, the 
relevant operator is F^. 

The unusual transformation law of the symbol reflects the fact that (1701) is not an auto- 
morphism. In the spectral theory of dynamical systems, one studies the dynamics of the 
billiard map (3 on B*Y through the associated 'Koopman' operator 

U : L\B*Y, da) L\B*Y, da), Uf{Q = f{m)- 

Here, da denotes the usual /5-invariant symplectic volume measure on B*Y . From the 
symplectic invariance it follows that W is a unitary operator. When (3 is ergodic, the unique 
invariant L^-normalized eigenfunction is a constant c. 

However, Egorov's theorem flTTl) in the boundary reduction involves the positive function 
7 G C°°(-B*y), and the relevant Koopman operator is 

Then T is not unitary on L'^{B*Y, da). Instead one has: 

• (i) The unique positive T- invariant function is given by 7. 
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• (ii) T is unitary relative to the inner product ((, )) on B*Y defined by the measure 

• (iii) When f3 is ergodic, the orthogonal projection P onto the invariant L^-eigenvectors 
has the form 



where c is as in flll2p . 

The proof of Theorem 16 . 21 then runs along similar lines to that of Theorem 14 . 1 1 but adjusted 
to the fact that that T is not unitary. 



This section is an exposition of Hassell's scarring result for the Bunimovich stadium. We 
follow [Has] and |Z7j . 

A stadium is a domain X = RUW gM.'^ which is formed by a rectangle R = [—a, a]x x 
[—P, P]y and where W = W-p U Wp are half-discs of radius (3 attached at either end. We fix 
the height /5 = 7r/2 and let a = tjd with t G [1,2]. The resulting stadium is denoted Xt. 

It has long been suspected that there exist exceptional sequences of eigenfunctions of X 
which have a singular concentration on the set of "bouncing ball" orbits of R. These are the 
vertical orbits in the central rectangle that repeatedly bounce orthogonally against the flat 
part of the boundary. The unit tangent vectors to the orbits define an invariant Lagrangian 
submanifold with boundary A C S*X. It is easy to construct approximate eigenfunctions 
which concentrate microlocally on this Lagrangian submanifold. Namely, let x{^) be a 
smooth cutoff supported in the central rectangle and form Vn = x{^) sinny. Then for any 
pseudo-differential operator A properly supported in X, 



J A 

where du is the unique normalized invariant measure on A. 

Numerical studies suggested that there also existed genuine eigenfunctions with the same 
limit. Recently, A. Hassell has proved this to be correct for almost all stadia. 

Theorem 7.1. The Laplacian on Xt is not QUE for almost every t G [1,2]. 

We now sketch the proof and develop related ideas on quantum unique ergodicty. The 
main idea is that the existence of the scarring bouncing ball quasi-modes implies that either 

• There exist actual modes with a similar scarring property, or 

• The spectrum has exceptional clustering around the bouncing ball quasi-eigenvalues 

Hassell then proves that the second alternative cannot occur for most stadia. We now 
explain the ideas in more detail. 

We first recall that a quasi- mode {ipk} is a sequence of L^-normalized functions which 
solve 




7. Hassell's scarring result for stadia 




|(A-/i^)^,||^. = 0(1), 
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for a sequence of quasi-eigenvalues /x^. By the spectral theorem it follows that there must 
exist true eigenvalues in the interval — K,fil + K] for some K > 0. Moreover, if E^^k 
denotes the spectral projection for A corresponding to this interval, then 

\\Ek,K^k - ^fc||L2 = 0{K~^). 

To maintain consistency with ([6]), i.e. with our use of frequencies fik rather than energies fij., 
we re-phrase this in terms of the projection E^^k for \/A in the interval ^/ jj^l + K]. 

For fixed K, this latter interval has width — . 

Given a quasimode {ipk}, the question arises of how many true eigenfunctions it takes to 
build the quasi-mode up to a small error. 

Definition: We say that a quasimode {ipk} of order with llV'fcllia = 1 has n{k) essential 
frequencies if 

n(fc) 

(72) ipk = ^CkjiPj + rik, ||?7fc||i2 = o(l). 

To be a quasi-mode of order zero, the frequencies Aj of the ipj must come from an interval 

[''^'^ ~ '^'^ ~'~ Mfc]' Hence the number of essential frequencies is bounded above by the 
number n{k) < N{k, y) of eigenvalues in the interval. Weyl's law for a/A allows considerable 
clustering and only gives N{k, = o{k) in the case where periodic orbits have measure zero. 
For instance, the quasi-eigenvalue might be a true eigenvalue with multiplicity saturating 
the Weyl bound. But a typical interval has a uniformly bounded number of A-eigenvalues 
in dimension 2 or equivalently a frequency interval of with 0{-^) has a uniformly bounded 
number of frequencies. The dichotomy above reflects the dichotomy as to whether exceptional 
clustering of eigenvalues occurs around the quasi-eigenvalues of A or whether there is a 
uniform bound on N{k, S). 

Proposition 7.2. // there exists a quasi-mode {ipk} of order for A with the properties: 

• (%) n{k) < C, V k; 

• (ii) {Atpk,ipk) js'M^Ad^i where dfi ^ dfi^. 
Then A is not QUE. 

The proof is based on the following lemma pertaining to near off-diagonal Wigner distri- 
butions. It gives an "everywhere" version of the off-diagonal part of Theorem 14.11 (2). 

Lemma 7.3. Suppose that is ergodic and A is QUE. Suppose that {(Aj,,, Xj^), ir 7^ jr} is a 
sequence of pairs of eigenvalues of ^/A such that Xi^ ~Xj^ ^ as r ^ oo. Then dWi^j^ 0. 

Proof. Let {Aj, Aj} be any sequence of pairs with the gap Aj — Xj ^ 0. Then by Egorov's 
theorem, any weak* limit du of the sequence {dWij} is a measure invariant under the geodesic 
flow. The weak limit is defined by the property that 

(73) {A*A^i,^,) ^ / WAl^du. 

Js*M 

If the eigenfunctions are real, then du is a real (signed) measure. 
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We now observe that any such weak* hmit must be a constant multiple of Liouville measure 
dfiL. Indeed, we first have: 

(74) \{A*A^,,y,,)\ < \{A*A^,,^,)\'/' 

Taking the limit along the sequence of pairs, we obtain 



(75) I / \aA\ di^\ < / \aA\ d^L- 

Js*M JS'M 

It follows that du « d^L (absolutely continuous). But dfii is an ergodic measure, so if 
du = fdfiL is an invariant measure with / G L}{d^i), then / is constant. Thus, 

(76) du = CdfiL, for some constant C. 

We now observe that C = if fi-Lipj (i.e. if i 7^ j). This follows if we substitute A = I in 
(1731) . use orthogonality and (1761) . 

□ 

We now complete the proof of the Proposition by arguing by contradiction. The frequencies 
must come from a shrinking frequency interval, so the hypothesis of the Proposition is 
satisfied. If A were QUE, we would have (in the notation of fl72l) : 



{A^k, ^k) = Zl"j=i CkjCkiiAipi, Lpj) + 0(1) 

= EJ=i 4j{^'^j^ Vj) + ESj-Li CkjCki{Aipi, ifj) + 0(1) 

= jg,j^j(TAdHL + 0{l), 

by Proposition 17.31 This contradicts (ii). In the last line we used Eji^i kfcjP = 1 + o(l), 
since HV^fcllia = 1. 

QED 



7.1. Proof of Hassell's scarring result. We apply and develop this reasoning in the case 
of the stadium. The quasi-eigenvalues of the Bunimovich stadium corresponding to bouncing 
ball quasi-modes are independently of the diameter t of the inner rectangle. 

By the above, it suffices to show that that there exists a sequence Uj 00 and a constant 
M (independent of j) so that there exist < M eigenvalues of A in [ra^ — K,v?j + K]. An 
somewhat different argument is given in |Hasj in this case: For each rij there exists a nor- 
malized eigenfunction u^. so that {ukpVk^) > \J \M. It suffices to choose the eigenfunction 

with eigenvalue in the interval with the largest component in the direction of Vky There 
exists one since 

3 

\\E[n^^K,n'^+KVn\ \ > ^• 

The sequence {un^,} cannot be Liouville distributed. Indeed, for any e > 0, let A be a 
self-adjoint semi-classical pseudo-differential operator properly supported in the rectangle so 
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that o"A < 1 and so that \\{Id — A)vn\ \ < e for large enough n. Then 

I 1 2 



Choose a sequence of operators A such that — ^)fn|| and so that the support of 
aA shrinks to the set of bouncing ball covectors. Then the mass of any quantum limit of 
{un^} must have mass > |M on A. 

Thus, the main point is to eliminate the possibility of exceptional clustering of eigenvalues 
around the quasi-eigenvalues. In fact, no reason is known why no exceptional clustering 
should occur. Hassell's idea is that it can however only occur for a measure zero set of 
diameters of the inner rectangle. The proof is based on Hadamard's variational formula for 
the variation of Dirichlet or Neumann eigenvalues under a variation of a domain. In the 
case at hand, the stadium is varied by horizontally (but not vertically) expanding the inner 
rectangle. In the simplest case of Dirichlet boundary conditions, the eigenvalues are forced 
to decrease as the rectangle is expanded. The QUE hypothesis forces them to decrease at 
a uniform rate. But then they can only rarely cluster at the fixed quasi-eigenvalues n^. If 
this ever happened, the cluster would move left of and there would not be time for a new 
cluster to arrive. 

Here is a more detailed sketch. Under the variation of Xt with infinitesimal variation 
vector field pt, Hadamard's variational formula gives, 

dEj{t) _ 



dt 

Then 



Pt{s){drrUj{t){s)) ds. 



dXt 



"f^ JdXt 



'dXt 

Let A{t) be the area of St- By Weyl's law, Ej{t) ~ ^~a(£)- Since the area of Xt grows 

linearly, we have on average Ej C-^y Theorem O gives the asymptotics individually 

for almost all eigenvalues. Let 

/.W= ! Pt{s)\u){t-s)\'ds. 



JdXt 

Then Ej = —Ejfj. Then Theorem 16.21 implies that — * weakly on the boundary 
along a subsequence of density one. QUE is the hypothesis that this occurs for the entire 
sequence, i.e. 

k f 



Hence, 



dSt 



^ = -kAit)il + o{l)), j^oo. 



Hence there is a lower bound to the velocity with which eigenvalues decrease as A(t) increases. 
Eigenvalues can therefore not concentrate in the fixed quasi-mode intervals [ri^ — K,n'^ + K] 
for all t. But then there are only a bounded number of eigenvalues in this interval; so 
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Proposition 17.21 implies QUE for the other Xt. A more detailed analysis shows that QUE 
holds for almost all t. 



8. Matrix elements of Fourier integral operators 

Our main focus in this survey goal is on the limits of diagonal matrix elements ( |T2l) 
of pseudo differential operators A or order zero relative to an orthonormal basis {<fj} of 
eigenfunctions. Difficult as it is, the limiting behavior of such matrix elements is one of 
the more accessible properties of eigenfunctions. To obtain more information, it would be 
useful to expand the class of operators A for which one can study matrix elements. In ^ 
we were essentially expanding the class from classical polyhomogeneous pseudo-differential 
operators to those in the "small-scale" calculus, i.e. supported in tubes around closed 
geodesies. Another way to expand the class of A is to consider matrix elements of Fourier 
integral operators. Motivation to consider matrix elements of Fourier integral operators 
comes partly from the fact that Hecke operators are Fourier integral operators, and this 
section is a preparation for the next section on Hecke operators ^ Other examples arise 
in quantum ergodic restriction theorems. Details on the claims to follow will be given in a 
forthcoming paper. 

By an FIO (Fourier integral operator) we mean an operator whose Schwartz kernel of F 
may be locally represented as a finite sum of oscillatory integrals, 

for some homogeneous phase and amplitude a. The only example we have seen so far in 
this survey is the wave group and of course pseudo-differential operators. 
The phase is said to generate the canonical relation 

C = {{x, <^;, y, -^'y) : <^^(x, y, ^) = 0} C T*M x T*M. 

The oscillatory integral also determines the principal ap of F, a 1/2 density along C. The 
data (C, ap) determines F up to a compact operator. In the case of the wave group, C is 
the graph of and the symbol is the canonical volume half-density on the graph. In the 
case of 'i/'DO's, the canonical relation is the diagonal (i.e. the graph of the identity map) 
and the symbol is the one defined above. We denote by /"(M, xM, C) the class of Fourier 
integral operators of order zero and canonical relation C. 

Hecke operators are also FIO's of a simple kind. As discussed in the next section, in that 
case C is a local canonical graph, i.e. the projections 

TTx : C ^ T*M, TTy : C ^ T*M 

are branched covering maps. Equivalently, C is the graph of a finitely many to one cor- 
respondence X '■ T*M —>■ T*M. Simpler examples of the same type include (i) F = is 
translation by an isometry of a Riemannian manifold (M, g) possessing an isometry, or (ii) 
Tf[x) = f{,9jx) + f{gj^x) on S'" corresponding to a finite set {gi, . . . ,gk} oi isometrics. 
The latter has been studied by Lubotzky-Phillips-Sarnak. 
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8.1. Matrix elements as linear functionals. The matrix elements 

(77) p,{F):={F^,,^,) 
define continuous linear functionals 

(78) Pj- : J°(M X M,C) ^ C 

with respect to the operator norm. It is simple to see that any limit poo of the sequence of 
functionals pj are functionals only of the principal symbol data (C, ap) which is bounded by 
the supremum of ap- Thus, as in the pseudo-differential case, poo is a measure on the space 
Vt^l"^ of continuous half-densities on C. 

The limiting behavior of pj (F) depends greatly on whether [F, -\/A] = (or is of negative 
order) or whether it has the same order as F. In effect, this is the issue of whether the 
canonical relation C underlying F is invariant under the geodesic flow or not. For instance, 
if C is the graph of a canonical transformation x and if x nioves the energy surface S*M = 
{(x,^) : \^\g = 1} to a new surface xi^gM) disjoint from S*M then {Fipj,{pj) 0. On 

the other hand, {Ae'^^^ipj, (pj) = e^^^^ {Aipj, ipj) so one gets the same quantum limits for 
F = y4e**^ as for A but must pick sparser subsequences to get the limit. 

Hecke operators have the property that [F, A] = 0. This implies that the canonical relation 
C is invariant under In general this is the only condition we need to study weak limits. The 
eigenfunction linear functionals pk are invariant in the the sense that pk{U~^FU^) = pk{F). 

Proposition 8.1. Let C C T*M x T*M be a local canonical graph, equipped with the 
pull-back of the symplectic volume measure on T*M. Assume C is invariant under g^ . Let 
F G /°(M X M,C) and identify its symbol with a scalar function relative to the graph half- 
density. Let Poo be a weak* limit of the functionals pj on P{M x N,C). Then there exists 
a complex measure v on SC = C f\ S*M x S*M of mass < 1 such that 



Poo(F) = 




which satisfies g^dv = dv. 

The proof is similar to that in the ipJ^O case. 

We note that /°(M x M, C) is a right and left module over ^°(M). Given F G /°(M x 
M, C) we consider the operators AF.FA G J°(M x M, C) where A G \I'°(X). We obtain a 
useful improvement on Proposition 18. 1[ 

Proposition 8.2. With the same hypotheses as in Proposition \8.1\ Assume further that F 
is self-adjoint and that pj{AF) ~ pj{FA). Then we have: 

Poo(7r^ao") = pooiTTyaa). 

This holds since 

poo(vr^acr) ~ pj{AF) = {AFLpj,Lpj) = {F*A*(pj, (pj) = {FA*(pj,(pj) ~ nyaa 

since a a* ~ and since op = ap. 

Let us consider some simple examples. First, suppose that F = Tg + T* where g is an 
isometry. Then C is the union of the graph of the lift of g to T*M (by its derivative) 
and the inverse graph. The above Proposition then says that g^v = v. This is obvious by 
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Egorov's theorem, and we can regard the Proposition as a generahzation of Egorov's theorem 
to symplectic correspondences. As a second example, let F = U*". Then {AFipj,ipj) = 
e^'^^^ {Aipj,ipj), so we see that sequences with unique quantum limits are sparser than in the 
pseudo-differential case. 

9. QUE OF Hecke eigenfunctions 

Our next topic is Hecke eigenfunctions. A detailed treatment would have to involve adelic 
dynamics, higher rank measure rigidity in the presence of positive entropy [LINDt ILIND2] 
and L-functions |HSl ISoundlj . Both areas are outside the scope of this survey. Fortunately, 
Lindenstrauss has written an expository article on adelic dynamics and higher rank rigid- 
ity for an earlier Current Developments in Mathematics |LIND2j . and Soundararajan has 
recently lectured on his L-function results in a 2009 Clay lecture. Sarnak has recently writ- 
ten exposition of the new results of Soundararajan and Holowinsky-Soundararajan results is 

So we continue in the same vein of explaining the microlocal (phase space) features of the 
Hecke eigenvalue problem. The distinguishing feature of Hecke eigenfunctions is that they 
have a special type of symmetry. As a result, their quantum limit measures have a special 
type of symmetry additional to geodesic flow invariance. The analysis of these additional 
symmetries is a key input into Lindenstrauss' QUE result. The main purpose of this section is 
to derive the exact symmetry. In fact, it appears that the exact symmetry was not previously 
determined; only a quasi-invariance condition of |RS] has been employed. 

We begin by recalling the definition of a Hecke operator. Let F be a co-compact or cofinite 
discrete subgroup of G = PSL{2, M) and let Xp be the corresponding compact (or finite area) 
hyperbolic surface. An element g E G,g is said to be in the commensurator Comm^T) 
if 

T'ig) :=Tng-'Tg 
is of finite index in F and g^^Tg. More precisely, 

d 

r = U F'(^)7j, (disjoint), 

i=i 

or equivalently 

d 

TgT = l^Tttj, where aj = gjj. 
i=i 

It is also possible to choose aj so that TgT = ljj=i '^i^. We then have a diagram of finite 
(non-Galois) covers: 

r'(^7)\H 

(79) TT / \p 

F\H ^ g-'rg\H. 

Here, 

'^i^'i9)z) = Tz, p{T'{g)z) = gTg~^{g-fjz) = gTz, 
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where in the definition of p any of the 7j could be used. The horizontal map is 2; — > g~'^z. 
A Hecke operator 

d 

(80) TJ{x) = figj.x) : L'{Xr) ^ ^^(Xr) 

i=i 

is the Radon transform p*7r* of the diagram. We note as a result that 

Tg-iu{z) = u{g~^z) 

takes F-invariant functions to (yf^^r^'-invariant functions. The Hecke operator is thus an 
averaging operator over orbits of the Hecke correspondence, which is the multi-valued holo- 
morphic map 

Cg{z) = {aiz, . . .,adz}. 

Its graph 

Tg = {{z,ajz) :zeXr)} 

is an algebraic curve in Xr x g~^Xg. 

The Hecke operators Tg form a commutative ring as g varies over Comm{V). By a Hecke 
eigenfunction is meant a joint eigenfunction of A and the ring of Hecke operators: 

TgUj = ^^j{g)uj, Auj = -X^jUj. 

9.1. Quantum limits on r'(g()\H. We now considered symmetries of quantum limits of 
Hecke eigenf unctions. We assume Uj is a Hecke eigenfunction. It is convenient to modify the 
definitions of the functionals pk as follows: 



{■K''Uj,p*Uj) {'K*Uj,p*Uj) ' 



(81) 



{A-K*Uj ,1X*Uj) 
{■K*Uj,p*Uj) ■ 

{An'^Uj ,p*Uj) 

{p*Uj,p*Uj) 



on A e *°(r'(^)\H). Thus, if Aq e *°(Xr) we have. 



(82) 



{AoP*'K*Uj, Uj) 



= {AoUj,Uj). 



{p*'K*Uj, Uj) 

Since Hecke operators are FlO's, limit measures of {ATgUj,Uj) and {TgAuj,Uj) live on the 
canonical relation underlying Tg, which is the lift of the graph of the Hecke correspondence 
to T*Xr. It is convenient to restate the result by forming a second diagram 



(83) 



r, C Xr X g-'Tg\U 

T^i y \ pi 

r\H ^ g~'rg\H. 
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Here, ni, pi are the natural projections. In the case of arithmetic groups such as SL{2,Z), 
Comm{r) is dense in G, i.e. there are many Hecke operators. The map 

i : T'{g)\U ^ r, 

defined by 

is a local diffeomorphic parametrization, and ttil = it, pit = p. 
The general Proposition 18.21 then specializes to: 

Proposition 9.1. If uj' is a weak* limit of p'j, and if uj is the weak* limit on Xr then 
-n^uj' = p^uj' = UJ. 

We may reformulate the proposition as follows: if we lift measures and functions to the 
universal cover H, then the data of the Hecke limit measures are a finite set of real signed 
measures {dvk} on the different copies of S*M.. The latter must project to the former under 
both projections, so we get the exact invariance property, 

d d 

(84) z/ = ^ z/„^. = ^ aj^Pa^ , 

i=i j=i 

for some complex measures Uj on the j sheet of the cover. The sum of the masses of the 
Vj must add up to one. This is simpler to see for Hecke operators on S"^. There we also 
have some number r of copies of S'^ and r isometrics Qj and define the Hecke operator (of 
Lubotzky-Phillips-Sarnak) as in (ISO]) . Then a quantum limit measure on 3*3"^ is the 
projection of r limit measures on the copies of 3*3"^ under both projections. 

Previously, only a quasi-invariance property of quantum limit measures of Hecke eigen- 
functions was proved by Rudnick-Sarnak |RSj . and used in |BH IW3j and elsewhere. In the 
notation above, it may be stated as 

(85) v{E) < dY,Wj*^a,{E)\. 

j 

As a simple application of fl84|) (which was also clear from the quasi-invariance condition), 
we sketch a proof that a quantum limit measure z/ of a sequence of Hecke eigenfunctions for 
r = 3L{2, Z) cannot be a periodic orbit measure p^. If it were, it would have to be ^ Uj as 
above. Each uj is an invariant signed measure so by ergodic decomposition, each must have 
the form Cj/i^ + tj where ^ c j = 1, where Tj are singular with respect to p^ and J2j = 0- 
But also the images of uj under the aj must have the same property. This is only possible 
if (yj{'~f) = 7 up to a translation by /3 G P whenever Uj ^ 0. At least one of the Cj ^ 0. At 
this point, one must actually consider the elements aj in the Hecke operator. In the case 
of P3L{2, Z) they are parabolic elements plus one elliptic element. But only a hyperbolic 
element can fix a geodesic (its axis). 

It would be interesting to see if one can determine how the z/^ are related to each other. 
Lindenstrauss's QUE result (in the next section) implies that that (Tg)*z/ = z/, which appears 
to say that the Ua are all the same. An apriori proof of this would simplify the proof of 
QUE. 
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9.2. QUE results of Lindenstrauss, Soundararajan (and Holowinsky-Soundararajan). 

We now briefly recall the QUE results on Hecke eigenfunctions. 

Theorem 9.2. (E. Lindenstrauss [LIND] ) Let Xr be a compact arithmetic hyperbolic sur- 
face. Then QUE holds for Hecke eigenfunctions, i.e. the entire sequence of Wigner distribu- 
tions of Hecke eigenfunctions tends to Liouville measure. 

In the non-compact finite area case of F = PSL{2, Z), Lindenstrauss proved: 

Theorem 9.3. (E. Lindenstrauss |LINDj ) Let Xp be the arithmetic hyperbolic surface de- 
fined by a congruence subgroup. Then any weak* limit of the sequence of Wigner distributions 
of Hecke eigenfunctions is a constant multiple of Liouville measure (the constant may depend 
on the sequence). 

Note that the quantum ergodicity theorem does not apply in this finite area case. In the 
arithmetic case there is a discrete spectrum corresponding to cuspidal eigenfunctions and a 
continuous spectrum corresponding to Eisenstein series. The following was proved in |Z8j : 

Theorem 9.4. Let Xp be the arithmetic hyperbolic surface defined by a congruence subgroup, 
and let {^Pj} be any orthonormal basis of cuspidal eigenfunctions. Then for any pseudo- 
differential operator of order zero and with compactly supported symbol, 

where Nc{X) ~ |Xr|A^ is the number of cuspidal eigenvalues < A. 

Recent work of Soundararajan that the constant equals one, i.e. there is no mass leakage 
at infinity. Hence the theorem is improved to 

Theorem 9.5. (E. Lindenstrauss [LIND] . Soundarajan [Soundlj ) Let Xr be the arithmetic 
hyperbolic surface defined by a congruence subgroup. Then any weak* limit of the sequence 
of Wigner distributions of Hecke eigenfunctions is is normalized Liouville measure. 

Holowinsky ^Holj and Holowinsky-Soundararajan |HSj have proven QUE results for holo- 
morphic forms (i.e. for the discrete series; see §1.15.41) . However, the methods are entirely 
different; they are based on the theory of L-functions and Poincare series. The Hecke prop- 
erty is exploited through multiplicativity of Fourier coefficients. 

9.3. Hecke QUE and multiplicities. A natural question (raised at the Clay talk of 
Soundararajan) is whether QUE holds for any orthonormal basis of eigenfunctions of an 
arithmetic quotient if it holds for the Hecke eigenbasis and if the discrete spectrum of A has 
bounded multiplicities. The proof that this is true is a simple modification of Proposition 
17.21 and we pause to sketch it now. 

Proposition 9.6. Suppose that one orthonormal basis of A is QUE, and that the eigenval- 
ues have uniformly bounded multiplicities. Then all orthonormal bases are QUE. 

Proof: As in Proposition 17. 21 consider sequences {(Aj^, ir 7^ jr} of pairs of eigenvalues 
of with Aj^ = Xj^. We then have that dWi^j,. — > 0. 

This has implications for all orthonormal basis as long as eigenvalue multiplicities are 
uniformly bounded. We run the argument regarding modes and quasi-modes but use the 



RECENT DEVELOPMENTS IN MATHEMATICAL QUANTUM CHAOS 45 

QUE orthonormal basis in the role of modes and any other orthonormal basis {ipk} in the 
role of quasi-modes. 

Lemma 9.7. If the eigenvalue multiplicities are uniformly bounded and there exists a sequence 
of eig en functions {ipk} 

• (i) n{k) < C, V k; 

• (ii) {Aipk,i^k) js'M^^dlJi where dfi 7^ dfiL- 

Then no orthonormal basis of eig enf unctions of A is QUE. 

Proof. We argue by contradiction. If {^Pj} were QUE, we would have (in the notation of 

{Atjjk, tpk) = YJlfli CkjCkiiAipi, Lpj) + 0(1) 

= Is*M^AdfiL + 0{l), 

by Proposition 19.61 This contradicts (ii). In the last line we used Yl^=i kfcjP = 1 + o(l), 
since llV'/cllia = 1- 

□ 

10. Variance estimates: Rate of quantum ergodicity and mixing 

A quantitative refinement of quantum ergodicity is to ask at what rate the sums in The- 
orem |4]T](1) tend to zero, i.e. to establish a rate of quantum ergodicity. More generally, we 
consider 'variances' of matrix elements. For diagonal matrix elements, we define: 

(86) Va{X) ^ E \{A'Pvf^)-^iA)\'- 

In the off-diagonal case one may view \ {Aipi, v^j)^ as analogous to \ {Aipj, ipj) —u^A)]"^. How- 
ever, the sums in f l57|) are double sums while those of fl86|) are single. One may also average 
over the shorter intervals [A, A -|- 1]. 

10.1. Quantum chaos conjectures. It is implicitly conjectured by Feingold-Peres in |FPj 
(11) that 



2ttp{E) ' 

where Ca(t) = J^^e~'^'^*{V^aA, crA)dt. In our notation. A-,- = h~^Ej (with h = Xj^ and 
Ej = A2) and p{E)dE ~ dN{X). 

On the basis of the analogy between | {A{pi, v^j) P and | {Aipj, ipj) —uj{A)\'^, it is conjectured 
in [FP] that 

Va{X) 



Ca-w(A)i{^) 



X'^-HoliSl) ■ 

The idea is that = -^{Vi i Vj) have the same matrix element asymptotics as eigenfunc- 
tions when Aj — Xj is sufficiently small. But then 2(y4(y9+, ipJ) = {Aipi, ipi) — {A(pj, ipj) when 
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A* = A. Since we are taking a difference, we may replace each matrix element by {A(fi, ipi) 
by {Aipi,ipi) — uj{A) (and also for ipj). The conjecture then assumes that {Aipi,ipi) — uj{A) 
has the same order of magnitude as {Aipi, ipi) — {Aipj, ipj). 

10.2. Rigorous results for Hecke eigenfunctions. Let Xr be the arithmetic modular 
surface, with F = SL{2, Z), and let T„, n ^ 1 denote the family of Hecke operators. In addi- 
tion to the Laplace eigenvalue problem fl28l) . the Maass- Hecke eigenforms are eigenfunctions 
of the Hecke operators, 

T„v2 = A^(n)v3. 
In the cusp, they have Fourier series expansions 

^(z) = ^/y^ p^{n)K„i2Tc\n\y)e{nx) 

where Kir is the i^-Bessel function and Pip{n) = A(^(n)p^(l). A special case of the variance 
sums ([HSD is the configuration space variance sum. 



A,<A 



Following the notation of |LS2j and |Zhj .we put dpn = ¥^r„d.V. 
Let ip G C^o('S'*Xr) and consider the distribution of 



Ratner's central limit theorem for geodesic flows implies that this random variable tends 
to a Gaussian with mean and the variance given by the non-negative Hermitian form on 

' ' 



J-oo ir\sL(2,R) V V e 2 yy 

Restrict Vc to C^q(X). Different irreducible representations are orthogonal under the qua- 
dratic form. Maass cusp forms (p are eigenvectors of Vc-. If the Laplace eigenvalue is | + r^, 
then the Vc eigenvalue is 

Ipj" 1 ir \ 14 



^4 2 



27r|F(i -2r)|2 



In |LSl ILS2] , Luo and Sarnak studied the quantum variance for holomorhic Hecke eigen- 
forms, i.e. holomorphic cusp forms in Sk{T) of even integral weight k for F. In this setting, 
the weight of the cusp form plays the role of the Laplace eigenvalue. They proved that for 
(p,ip E C^q(X), as the weight K oo, 

E E ^(l,SymV)/i/(^)MVO 

where B{ip,il)) is a non-negative Hermitian form on C^q{X). B is diagonalized by the 
orthonormal basis of Maass-Hecke cusp forms and the eigenvalues of B at ifjj is '^LiJ^^ipj). 
For the notation L(l, Sym^) we refer to |LS2t IZhj . 
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In |Zh] . P. Zhao generalized the result in |LSt ILS2] to Maass cusp forms. Let <fj{z) be the 
j-th Maass- Hecke eigenform, with the Laplacian eigenvalues \ + r'^j- 

Theorem 10.1. [Zh] Fix any e > 0. Then we have 

(^g-l^j ^g-i^J )L(l,sym ¥?j)/ij-(¥?)7I7(?/') 

is diagonalized by the orthonormal basis {<fj} of Maass-Hecke cusp forms and the eigen- 
value ofV at a Maass-Hecke cusp form ip is 

It would be desirable to remove the weights. The unweighted version should be (P. Zhao, 
personal communication; to appear) 



(88) 5^/i,-(^)7I^.(^)~A\/(^,V^). 

A,<A 

10.3. General case. The only rigorous result to date which is valid on general Riemannian 
manifolds with hyperbolic geodesic flow is the logarithmic decay: 

Theorem 10.2. [Z4] For any {M,g) with hyperbolic geodesic flow, 



It is again based on Ratner's central limit theorem for the geodesic flow. The logarithm 
reflects the exponential blow up in time of remainder estimates for traces involving the wave 
group associated to hyperbolic flows and thus the necessity of keeping the time less than the 
Ehrenfest time (1181) . It would be surprising if the logarithmic decay is sharp for Laplacians. 
However, R. Schubert shows in [Schuj ISchu2] that the estimate is sharp in the case of two- 
dimensional hyperbolic quantum cat maps. Hence the estimate cannot be improved by 
semi-classical arguments that hold in both settings. 

10.4. Variance and Patterson-Sullivan distributions. When /g.Xp (^^hl = 0, the Luo- 
Sarnak variance results for Hecke eigenfunctions have the form, 

(89) ^ \{a,dW,r,)\'r^\via,a). 

Even if one knew that asymptotics of this kind should occur, it does not seem apriori obvious 
that the bilinear form on the right side should be invariant under g^ because V{a, a) occurs 
in the j term in the variance asymptotics and only only knows that the Wigner distributions 
are invariant modulo terms of order A quick way to see that the bilinear form must be 
(yf*-invariant is via the Patterson-Sullivan distributions. From ( 15T1) . we have 

(90) Mat 5Z K«'^^-.)l' = MAy 5Z \{a,clPS.r,)\' + 0{X-'). 
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Indeed, the difference is of the form Sj |rj|<A T'^^ i^y dPSir-) ■ dPSir^) + 0(A ^). 
The sum in the last expression is also of order 0(A~^) (e.g. one can go back and express each 
factor of PS with one of Wir- and apply the Luo-Sarnak asymptotics again). Since PS'j^ is 
(yf*-invariant, it follows that V{a, a) must be invariant. This argument is equally valid on 
any compact hyperbolic surface, but of course there is no proof of such asymptotics except in 
the arithmetic Hecke case. It would be interesting to draw further relations between V"(a, a) 
and PSiry 

11. Entropy of quantum limits on manifolds with Anosov geodesic flow 

So far, the results on quantum limits have basically used the symmetry (or invariance) 
properties of the limits. But generic chaotic systems have no symmetries. What is there to 
constrain the huge number of possible limits? 

The recent results of Anantharaman [A], Anantharaman- Nonnenmacher |ANj and Riviere 
|Riv] ■ give a very interesting answer to this question (see also |ANKj ) for (M, g) with Anosov 
geodesic flow. They use the quantum mechanics to prove lower bounds on entropies of 
quantum limit measures. The lower bounds eliminate many of the possible limits, e.g. they 
disqualify finite sums of periodic orbit measures. 

The purpose of this section is to present the results of and [AN] in some detail. Both 
articles are based on a hyperbolic dispersive estimate (called the main estimate in jA] ) 
which, roughly speaking, measures the norm of quantized cylinder set operators in terms 
of h and the length of the cylinder. But they use the estimate in quite different ways. In 
[A], it is used in combination with an analysis of certain special covers of S*M by cylinder 
sets that are adapted to the eigenfunctions. One of the important results is an estimate 
on the topological entropy /itop(supp(z/o)) of the support of any quantum limit measure. In 
[AN] the key tool is the "entropic uncertainty principle". It leads to a lower bound for the 
Kolomogorv-Sinai entropy hKsi^o) of the quantum limit. 

There now exist several excellent and authoritative expositions of the KS entropy bounds 
and the entropic uncertainty principle |A2t lANKt IAN2t ICV3j in addition to the well written 
initial articles jAllANj (which take considerable care to give intuitive explanations of technical 



steps). For this reason, we emphasize the approach in [A]. We closely follow the original 
references in discussing heuristic reasons for the lower bounds and outlining rigorous proofs. 
We also discuss earlier entropy lower bounds of Bourgain-Lindenstrauss |BLj and Wolpert 
|W3j in the case of Hecke eigenfunctions. 
Before stating the results, we review the various notions of entropy. 

11.1. KS entropy. We first recall the definition of the KS entropy of an invariant prob- 
ability measure fi for the geodesic flow. Roughly speaking, the entropy measures the av- 
erage complexity of yU-typical orbits. In the Kolmogorov- Sinai entropy, one starts with a 
partition V = {Ei, . . . ,Ek) of S*M and defines the Shannon entropy of the partition by 

hp{fj,) = Yl^=i f^i^j) log/^(-^j)- Under iterates of the time one map g of the geodesic flow, 
one refines the partition from the sets Ej to the cylinder sets of length n: 

:= {V2 := E^, n g-'E^, n ■ ■ ■ H g-'^+'E^,^,, : [«o, • • • , «n-i] e N"}. 

One defines hn{V, fi) to be the Shannon entropy of this partition and then defines hxsif^, — 
lim„^oo ^hn{n,V). Then hKsil^) = sup^ /iii-s(/i, P). 
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The measure /i(-Eao H 9~^Ea^ fl ■ ■ ■ fl g^'^^^Ea„_-^) is the probabihty with respect to /i that 
an orbit successively visits , Ea^ , ' ■ The entropy measures the exponential decay rate 
of these probabilities for large times. 

11.2. Symbolic coding and cylinder sets. In the {M,g) setting, we fix a partition {Mk} 
of M and a corresponding partition T*Mk of T*M. Let Pa^ be the characteristic function 
of T*Mk. (Later it must to be smoothed out). Let S = {1, . . . , i}^ where i is the number of 
elements of the partition V'^\ To each tangent vector v E S*M, one can associated a unique 
element I{v) = (ak) G S so that g"'v G Pa^ for all n G Z. This gives a symbolic coding map 
I : S*M Ti. The time one map then conjugates under the coding map to the shift 
a{{aj)) = ((aj+i)) on admissible sequences, i.e. sequences in the image of the coding map. 
An invariant measure thus corresponds to a shift invariant measure /io on S. 

A cylinder set [ao, . . . , On-i] C S of length n is the subset of S formed by sequences with 
the given initial segment. The set of such cylinder sets of length n is denoted S„. Cylinder 
sets for S are not the same as cylinder sets for g^, which have the form Pa^ fl g^^Pa^ fl ■ ■ ■ fl 

The Ho measure of a S- cylinder set is by definition, 

/io([ao, • • • , aj) = MPao n g~'Pa, n • ■ ■ n ^""+'Pa„_i). 

11.3. Bowen balls. Cylinder sets are closely related to Bowen balls, i.e. balls in the Bowen 
metric. For any smooth map f : X —>■ X, the Bowen metric at time n is defined by 

di{x, y) = max d{fx, fy). 

0<i<n— 1 

Let Bfi^{x,r) be the r ball around x in this metric. 

In the continuous time case of the geodesic flow g^, one defines the metric 

(91) dT{p,p') = maxd{g\p),g\p')). 

te[OT] 

The Bowen ball BT{p,r) at time T is the r- ball in this metric. 

As with cylinder sets, if p,p' lie in BT{p,e) then their orbits are e-close for the interval 
[0,T]. This is not quite the same as running through the same elements of a partition but 
for large T the balls and cylinders are rather similar. This is because the geodesic flow g'^ 
stretches everything by a factor of e* in the unstable direction, and contracts everything 
by e~* in the stable direction; it preserves distances along the geodesic flow. In variable 
curvature and in higher dimensions, the 'cube' becomes a rectangular parallelopiped whose 
axes are determined by the Lyapunov exponents. 

In the case of geodesic flows of compact hyperbolic manifolds of constant curvature — 1, 
Bowen balls -Bt(P; are roughly of radius re~^ in the unstable direction, and r in the stable 
direction and geodesic flow directions. In the hyperbolic case of G/T with co-compact F, the 
tube (in the notation of §1.15.11) is 

(92) BT{p,e)) = a((— r, r))n_(r, r))ri+((— e~"^r, e~"^r)). 

To make the ball symmetric with respect to the stable and unstable directions is to make 
it symmetric with respect to time reversal. One uses the time interval [— T/2,T/2] instead 
of [0, T] and defines the new distance, 

(93) d'Ap,p) = ^^m^x^^^^d'{g\p),g\p')), 
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We then denote the r ball by B'rp[p, r). In the hyperbolic case, 

(94) B'Ap,e)) = a{{-e,e))n^{{-e^^''e,e-^'h))n^{{-e-^l\,e-^'h)). 

Thus, the Bowen ball in the symmetric case (in constant curvature) is a ball (or cube) of 
radius re~'^ in the transverse direction to the geodesic flow. The length along the flow is not 
important. 

11.4. Brin-Katok local entropy. Kolmogorov- Sinai entropy of an invariant measure is 
related to the local dimension of the measure on Bowen balls. This is the approach in 
Bourgain-Lindenstrauss |BL] and in [LIND] , and stems from work of Young and Ledrappier- 
Young. 

Define the local entropy on an invariant measure yU for a map / by 

KU.x) = limlimsup^^^^i^^^^^^^ = limliminf-^"^^^'^"^^'^^^ 

5^0 n->oo n <5^0 n^oo n 

Brin-Katok proved that both limits exist, that the local entropy hn{f,x) is / invariant 
and 

hM) = / h^{f,x)dfx{x). 
J X 

The definition and result is similar in the geodesic flow case. The local entropy of /i at p 
is defined by, 

limliminf-;^log/i(5T(p,e)) = limlim sup log /x(St(p, e)) := hKs{f^,p)- 

11.5. Ergodic decomposition. Let v be an invariant measure. It may be expressed as a 
convex combination of ergodic measures, which are extreme points of the compact convex 
set of invariant measures. There is a concrete formula, 

(95) u = j uldv{x) 

for its ergodic decomposition, where z/f is the orbital average through x G G/T. It is a fact 
that z/f ergodic for i/-a.e. x. Then the local entropy is 

a fi(T(x,e)) 

(96) /ia(t)(Px) = lim — i , a.s.p. 

e^O log e 

The Brin-Katok theorem states that the global KS entropy is 

h{^) = I h{ul)dv{x). 



11.6. Topological entropy of an invariant subset. Anantharaman's first entropy bound 
refers to the topological entropy of an invariant set F. Then, by definition, 

htop{F) < A <^=^ V(5 > 0, 3C > : F can be covered by at most 

(97) 

(j^n{\+S) cylinciers of length n, Vn. 
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11.7. Bourgain-Lindenstrauss entropy bound. Before discussing the work of Anantharman- 
Nonnenmacher, we review a simpler entropy bound of Bourgain-Lindenstrauss. After initial 
work of Wolpert |W3] , Bourgain-Lindenstrauss jBL] obtained a strong lower bound on en- 
tropies of quantum limit measures associated to Hecke eigenfunctions. Although it was later 
surpassed by Lindenstrauss' QUE theorem |LINDj . it is simple to state and illustrates the 
use of the local entropy and local dimensions to measure the KS entropy. 

The setting is a compact or finite area arithmetic hyperbolic quotient. We have not 
reviewed local entropy in the non-compact finite area case, but proceed by analogy with the 
compact case. Then 

B{e,To) == a((-ro,ro))n_(e,e))n+((-e,e)) 

is a Bowen ball around the identity element and xB{e, tq) is its left translate hy x E G. Here 

e = re~'^l'^ . 

Theorem 11.1. jBL] Let Y C SL{2, Z) he a congruence subgroup. Let n be a quantum limit. 
Then for any compact subset K C T\G and any x E K , 

l,{xB{e,To))<e'/'. 

In the flow-box notation above, 
(98) /i„(5T(p,e)) <Ce-^/^ 

This implies that any ergodic component of any quantum limit has entropy > |. 

Corollary 11.2. [BL] Almost every component of a quantum limit measure has entropy 
> 1/9. The HausdorfJ dimenson of the support of the limit measure is > 1 -f 1/9. 

Theorem 111.11 is a consequence of the following uniform upper bound for masses in small 
tubes around geodesic segments in configuration space: 

Theorem 11.3. |BLj Let T C SL{2, Z) be a congruence subgroup. Let fi be a quantum limit. 
Then for any compact subset K C r\G and any x E K, and for any Hecke eigenfunction 

JxB{r,To) 

One may rewrite this result in terms of a Riesz-energy of the quantum limit measures: 

Corollary 11.4. Let ^ be a quantum limit of a sequence of Hecke-Maass eigenfunctions. 
Then for n < 2/9, 

f r dfi{x)dfi{y) 
JmJm dM{,x,yY+^ 

Although these results were superseded by Lindenstrauss' QUE bound, we briefly sketch 
the idea of the proof since it makes an interesting contrast to the Anantharaman-Nonnenmacher 
bound. They prove that there exists a set W of integers of size e~^/^ so that for n E W one 
has (roughly) that the translates yB{e, tq) of the small balls by y in the Hecke correspondence 
for Tn are all pairwise disjoint. Since 
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one has 

= E.eT„(x)/.B(e,ro)l<^A(2;)|W. 

Now sum over n E W and use the disjointness of the small balls zB{e, tq) and the 
normalization of the eigenfunction to obtain 

11.8. Anantharman and Anantharaman-Nonnenacher lower bound. We now turn 
to the lower bounds on entropy in the general Anosov case given in [X] and in [ANt lANKj . 
The mechanisms proving the lower bounds and the results are somewhat different between 
the two articles. The result of [A] gives a lower bound for the topological entropy htopisuppuo) 
of the support of a quantum limit while the main result of [ANj gives a lower bound on the 
metric or Kolmogorov Sinai entropy hg = hxs of the quantum limits. For the definition of 
htop{F), see dnH). 

Theorem 11.5. [A] Let z/q = J^^^j UqcIuq^x) be the ergodic decomposition ofv^. Then there 
exists K > and two continuous decreasing functions r : [0, 1] [0, 1] and 'd : (0, 1] M+ 
with r(0) = 1,'(9(0) = oo, such that 



^o[{x : Kiyl) > ^(1 - 5) j] > (-^)2(1 - r{6)). 

Hence 

• hg{uo) > 0; 

• htopisupp Z/q) > f . 

It follows that a positive proportion of the ergodic components of z/q must have KS entropy 
arbitrarily close to ^. 

This is proved as a consequence of a Proposition which may have other applications: 

Theorem 11.6. [A] (Proposition 2.0.4) With the same notation as in Theorem \11.5\ and 
Remark ] 1 1 . 1 (A Let F be a subset with htop{F) < ^{1 — 6). Then, 

MF)<[l-{^^ryi-r{6))<l. 

Hence supp uq cannot equal F. 

In [A] Theorem 1.1.2, there is a generalization to quasi-modes which introduces a constant 
c. For simplicity we only consider eigenfunctions. 

In the subsequent article of Anantharaman-Nonnenmacher, [ANt lANK] the authors obtain 
a quantitative lower bound on the KS entropy: 

Theorem 11.7. |ANl lANK] Let ii be a semiclassical measure associated to the eigenfunc- 
tions of the Laplacian on M . Then its metric entropy satisfies 

id-1) 



(99) hKsif^) > 



logJ"(p)d/i(p) 

S*M 



"A max ; 
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where d = dimM. Amax = liHi|f|^oo 7 logsup^g^- Idg^l is the maximal expansion rate of the 
geodesic flow on S (l^j and J" is the unstable Jacobian determinant 

When M has constant sectional curvature —1, the theorem imphes that 

(100) hKsif^) > 
They state the conjecture 

Conjecture 11.8. |AN] Let {M,g) have Anosov geodesic flow. Then for any quantum 
limit measure uq, 

[ logriv)duo{v) . 

J S*M 

This conjecture has recently been proved by G. Riviere for surfaces of negative curvature 
[Rk]. 

This conjecture does not imply QUE. In fact the counterexamples in |FNBt IFN] to QUE 
in the case of "quantum cat maps" satisfy the condition of the conjecture. So there is some 
chance that the conjecture is best possible for (M, g) with ergodic geodesic flow. 

To the author's knowledge, the corollary that quantum limits in the Anosov case cannot 
be pure periodic orbit measures has not been proved a simpler way that by applying the 
theorems above. In §11.211 we will go over the proof in that special case and see how it ties 
together with mass concentration of eigenf unctions around hyperbolic closed geodesies in §5] 

11.9. Problem with semi-classical estimates and Bowen balls. As a first attempt to 
estimate KS entropies of quantum limits, one might try to study the local entropy formula 
as in Bourgain-Lindenstrauss. Let -Bt(p, e) be the small tube around p G S'*M, i.e. of length 
e in the flow direction and e~"^/^ in the stable/unstable directions. As in the Bourgain- 
Lindenstrauss bound, one would like to understand the decay of fin{BT{p, e)) as T 00 and 
n — >■ 00. In the brief expository article |AN2] . the authors give the heuristic estimate 

(101) /i„(5T(p,e)) < CAf^e"^^, d = dimM 

on the local dimensions of the quantum measures for a sequence of eigenfunctions. Unlike 
the uniform Bourgain-Lindenstrauss estimate, the estimate is A„-dependent and is only non- 
trivial for times T > logA„. But then the ball radius (in the stable/unstable directions) is 
= An^^^, which is just below the minimal scale allowed by the uncertainty principle ( §1.4p . 
To obtain a useful entropy bound one needs to take T = 21ogA„ and then e"^ = A~^, well 
below the Heisenberg uncertainty threshold. So although it is attractive, the local entropy 
is difficult to use. Recall that Bourgain-Lindenstrauss used the many elements of the Hecke 
correspondences to move the small ball around so that it almost fills out the manifold, and 
then used the normalization to estimate the total mass. There doesn't seem to exist a 
similar mechanism to build up a big mass from the small balls in the general Anosov case. 

11.10. Some important parameters. Two important parameters K,,'d appear here and 
in the statements and proofs in [A]. They represent time scales relative to the Ehrenfest 
time: 

• k: Semi-classical estimates only work when n < n\ log^|, i.e. k\ logh\ plays the role 
of the Ehrenfest time. 
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• t?: The main estimate is only useful when h ^I'^e "^/^ << 1 or n > -(91 log/i|. 

11.11. Cylinder set operators in quantum mechanics. To get out of the too-small 
Bowen ball impasse, Anantharaman [X] and Ananthraman-Nonnenmacher study "quantum 
measures" of cylinder sets C rather than small balls. The emphasis (and results) on quan- 
tum cylinder sets is one of the key innovations in The classical cylinder sets are not 
directly involved in the main estimates of [XI lANj . Rather they are quantized as cylinder set 
operators. Then their "quantum measure" or quantum entropy is studied. To distinguish 
notationally between classical and quantum objects, we put a ■ on the quantum operator. 

To define quantum cylinder set operators, one quantizes the partition to define a smooth 
quantum partition of unity Pk by smoothing out the characteristic functions of M^. Let 



a = [ao, . . . , an-i G S„ be a cylinder set in sequence space S. The corresponding quantum 
cylinder operator is 

(102) Pa = Pa„.An - l)Pa„-2(n - 2) ■ ■ ■ P„„ 
where a = [ao, ■ ■ ■ , On-i] where 

(103) P{k) = U*^PU^. 

Here, U = is the propagator at unit time (or in the semi-classical framework, U = e*^"^/^ 



(see [A] (1.3.4)). 

The analogue of the measure of a cylinder set in quantum mechanics is given by the 
matrix element of the cylinder set operator in the energy state, transported to S. We use 
the semi-classical notation of [X], but it is easily converted to homogeneous notation. 

Definition: (See jS], (1.3.4)) Let ipn be an eigenfunction of A. Define the associated 
quantum "measure" of cylinder sets C = [ao, . . . , a„_i] G S„ by 

(104) /i?i([ao, • • •,"«]) = {PaM ■ ■ ■ PaQ{0)i'h,i'h)- 
Thus, one "quantizes" the cylinder set C = [ao, ■ ■ ■ , On-i] as the operator 

(105) C = f/-("-^)p„_,f/P„„„, ■ ■ ■ UPa, 
For each eigenfunction, one obtains a linear functional 

(106) fXFi{[ao, . . . , an^i]) = pn{C) 

in the notation of states of ^ The 'quantum measures' are shift invariant, i.e. 

Hh{[ao, ■ ■ ■ ,an-i]) = Hh{<y~^[ao, ■ ■ . ,an-i])- 

This says, 

which is true since Ph{A) = {A(pfi,(ph) is an invariant state. 

The matrix element {Paipn, '4^h) is the probability (amplitude) that the particle in state ipn 
visits the phase space regions Pq,q, Pa^, . . . , Pa^-i at times 0, 1, . . . , n — 1. 

Thus, the states pn or alternatively the Wigner distributions of ipn are transported to S 
to define linear functionals yU/j on the span of the cylinder functions on E. They exactly 
invariant under the shift map. 
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The functionals /x^ are not positive measures since Cn is not a positive operator. However 
as in the proof of the quantum ergodicity Theorem 14.11 one has the Schwartz inequahty, 

(107) \pn{C)\' < pniCX), 
reflecting the positivity of the state pn- 

11.12. Main estimate for cylinder sets of S. The so-called main estimate of [2] implies 

that for every cylinder set C G S„, one has 

(108) \pn{C)\ < C^/i-'^/2g-nA/2^^ ^ 0{e))'', uniformly for n < (3 \ logh\. 

Here, (3 can he taken arbitrarily large. This estimate is similar in spirit to (110 II) but cir- 
cumvents the uncertainty principle. It is one of the essential quantum mechanical (or semi- 
classical) ingredient for getting a lower bound on htopisuppuo) of supports of quantum limit 
measures (see §11.151) . 

11.13. Quantization of cylinder sets versus quantized cylinder operators. Let us 

explain how fllOSp gets around the uncertainty principle. A fundamental feature of quanti- 
zation is that it is not a homomorphism. Thus, the quantization of the (smoothed) char- 
acteristic function of the cylinder set Pao H g^^Pa^ fl • • ■ fl g^"'Pa„ is very different from 
Pa„-iUPa„_2 ' ' 'UPao, the Ordering which occurs in the quantum cylinder operator Pa in 
(11021) . The problematic Bowen ball estimate (llOip would arise if one quantized the Bowen 
ball (or cylinder set) in the first sense and then took its matrix element. Quantizing in the 
other order makes the main estimate possible. 

11.14. Quantized cylinder operators and covering numbers. As motivation to study 
quantum measures of cylinder sets, Anantharaman observes that the main estimate suggests 
an obstruction to a quantum measure concentrating on a set F with htop{F) < |-. For each n 

such a set can be covered by Ce"'^^''^~^^ cylinder sets of length n. By the main estimate (llOSp . 
each of the cylinder sets has quantum measure \phiC) \ < (7e-'^/2iog?i.-nA/2 uniformly for n < 
K I log^l for any k. As yet, the combination does not disprove that suppz/Q = F, i.e. z/o(F) = 
1, since the the large factor of /i""^/^ is not cancelled in the resulting estimate Ph{F) < 
(jg^n{^{i-5)-d/2iogh-nA/2 ^ additional idea is needed to decouple the large time parameter 
n into independent time parameters. Roughly speaking, this is done by a sub-multiplicative 
property. In §11.20[ the ideas in the above paragraph are made more precise and effective. 

11.15. Upper semi-continuity, sub-multiplicativity and sub-additivity. Before get- 
ting into the details of the (sketch of) proof, it is useful to consider why there should exist 
a non-trivial lower bound on hKsi^o) for a quantum limit measure? What constrains a 
quantum limit measure to have higher entropy than some given invariant measure? 

First, the classical entropy of an invariant is upper semi-continuous with respect to weak 
convergence. This is in the right direction for proving lower entropy bounds for quantum 
limits. If one can obtain a lower bound on the 'entropies' of the quantum measures 
then the entropy of the limit can only jump up. Of course, the quantum measures are not 
g*' invariant measures and so upper semi-continuity does not literally apply. One needs to 
define a useful notion of entropy for the quantum measures and prove some version of upper 
semi-continuity for it. The USC property of the classical entropy suggests that there should 
exist some analogue on the quantum level. 
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In [X], entropies are studied via special covers by cylinder sets. For these, there is a 
sub-multiplicative estimate on the number of elements in the cover. 

In |AN] . a quite different idea is used: the authors employ a notion of quantum entropy 
due to Maassen-Uffink and prove a certain lower bound for the quantum entropy called the 
entropy uncertainty inequality. It is the origin of the lower bound for the quantum entropies 
and the KS entropy in Theorem 111.71 (see §11.171) . As a replacement for the USC property 
of the entropy, a certain sub-additivity property for the quantum entropy is proved (see 

mm- 

In the next three subsections, we go over the main ingredients in the proof. Then we give 
a quick sketch of the full proof. 

11.16. Main estimate in more detail. The main estimate in [S] (Theorem 1.3.3) was 
refined in |ANt lANK] . To state the results, we need some further notation. Let nE{h) denote 
the Ehrenfest time 

r (i-^)iiognr 

riEin) = 

^max 

where the brackets denote the integer part and 5 is a certain small number arising in energy 
localization (for the definition of Amax see (jS])). The authors also introduce a discrete 'coarse- 
grained' unstable Jacobian 

Jr («o,«i) := sup{ J"(p) : p G T*n^,nS' : g'p G T*fi,J, 

for ao, ai = 1, . . . , K. Here, Qj are small open neighborhoods of the partition sets Mj. For 
a sequence a = (ao, . . . , On-i) of symbols of length n, one defines 

J^(a) := Ji{ao, ai) ■ ■ ■ J{'(a„-2, a„_i). 

Theorem 11.9. (See Theorem 3.5 of |ANK] j Given a partition V^^^ and 5, 5' > small 
enough, there exists ^(o) ,5 5/ such that, for any h < ^(0)^5^5/, for any positive integer n < 
nEih), and any pair of sequences a, a' of length n. 



(109) |P:,[/"P,Op(x(")) I <C/l-(^-i+^^) VJ«(a)J«(«')- 
Here, d = dimM and the constants c, C only depend on {M,g). 

There is a more refined version using a sharper energy cutoff in [AN] . 

Theorem 11.10. Given a partition V^^\ k > and 5 > small enough, there exists 
f^v<-o),5,K such that, for any h < hp(o)^s,K> f^'^ ^'^y positive integer n < hi\ \ogh\, and any pair 
of sequences a, a' of length n, 

(110) |P*, f/" p„ Op(x^"^) I < c h-^^+^^ e-(^-^)"(i + 0{h^))'' . 

Here, d = dimM. 

To prove this, one shows that any state of the form Op(x^*'')^ can be decomposed as a 

(d-l) 

superposition of essentially h 2 normalized Lagrangian states, supported on Lagrangian 
manifolds transverse to the stable leaves of the flow. The action of the operator Pq, on such 
Lagrangian states is intuitively as follows: each application of U stretches the Lagrangian 
in the unstable direction (the rate of elongation being described by the unstable Jacobian) 
whereas each multiplication by Pa cuts off the small piece of the Lagrangian in the ath cell. 
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This iteration of stretching and cutting accounts for the exponential decay. A somewhat 
more detailed exposition is in §11.221 

This estimate can be reformulated as a statement about matrix elements of the operators 
Pa relative to eigenfunctions (see Theorem 3.1 of [ANK] ) 

Theorem 11.11. Let i\ = 1^- For any n > there exists > so that, uniformly for 
h < and all n < K \ log h\ , and for all a = (ao, . . . , ctn-i) ^ [1, k]" fl Z", 

\\PaM < 2(27r;i)^'^/V"^(l + 0(e))". 

The cutoff operator Op{x^"'^) is supported near the energy surface A^, and the eigenfunc- 
tions are supported on their energy surfaces. Here, A is the "smallest" expansion rate (see 

(H). 

We recall that {Pa'iph, "iph) is the probability (amplitude) that the particle in state ipn visits 
the phase space regions Pq,„, Pq,^ , . . . , Pa„_i at times 0, 1, . . . , — 1. The main estimates show 
that this probability decays exponentially fast with n, at rate ^. However, the decay only 
starts near the Eherenfest time rii = ^^ii^-^. 

11.17. Entropy uncertainty principle. The source of the lower bound for hKsi^o) 
is most simply explained in |ANj . It is based on the 'entropic uncertainty principle' of 
Maassen- Uffink. There are several notions of quantum or non- commutative entropy, but for 
applications to eigenfunctions it is important to find one with good semi-classical properties. 

Let {H, (., .)) be a complex Hilbert space, and let WipW = ■>/ ip) denote the associated 
norm. The quantum notion of partition is a family vr = {T^k)k=i,...,M of operators on Ti such 
that X^ILi ^fc^fc = ^d. If \ = 1, the entropy of with respect to the partition vr is define 
by 

k=l 

We note that the quantum analogue of an invariant probability measure fi is an invariant 
state p, and the direct analogue of the entropy of the partition would be ^ p^TTkul) log pijikT^^)- 
If the state is p{A) = {Aip,ip) then p{7ikK) = IkfcV'iP- 

The dynamics is generated by a unitary operator U on Ti. We now state a simple version 
of the entropy uncertainty inequality of Maasen- Uffink. A more elaborate version in (XJ lANt 
lANK] gives a lower bound for a certain 'pressure'. 

Theorem 11.12. For any e > 0, for any normalized ETi, 

K{U^)+K{^) > -21ogc(W), 

where 

c{U) = sup Kcfe, Uej)\ 

is the supremum of all matrix elements in the orthonormal basis {cj}. In particular, hj^lilj) > 
— logc(W) if ip is an eigenfunction ofU. 

The intuitive idea is that if a unitary matrix has small entries, then each of its eigenvectors 
must have large Shannon entropy. 
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This uncertainty principle is applied to the entropies of the partitions defined above. For 
the propagator, we put lA = e*^^^ is the wave operator at the 'Ehrenfest time' = 
Or in the semi-classical framework (where h = j), one can use the Hamiltonian is H = h'^A 

and the time evolution U = e'"-'^^^^''^ with nE(h) = 
Applied to an eigenfunction of W, one has 

(111) KM= J2 \\P:A\\'^og\\p*M\^ 

\a\=nE 

and obviously 
On the right side, 

(112) c(W) = max | |P„/?7"^P„0p(x^"^^)| I 

\a\=\a'\=nE 

where x^"^'' is a very sharp energy cutoff supported in a tubular neighborhood £^ := H^^(l — 
e, 1 + e) of ^ = S*M of width 2h^-^e''^ for a given S > 0. 

11.18. Sub-additivity. Another key component in the proof is that the quantum entropy 
is almost sub- additive for ^ < | logt|. 

Sub-additivity of the classical KS entropy follows from the concavity of the function 
— xlogx. It states that the sequence {/i„(/i,P)} is sub-additive, i.e. 

hn+m {f^,P) < hn{fi,P) + h^{fi,P). 

However, this is false for the quantum measures. 

The correct statement is as follows: There exists a function R{nQ, K) such that lim^^^o Ri^o, K) 
for all uq ^ 7j and such that for all rio, "n, G N with riQ + n = ilz|)ii2ii!l ^ and for any normalized 

'^max 

eigenfunction (fx, one has 

(113) hp^o+^ifx) < h-pnoi^ipx) + h'pn{ipx) + R{no,h). 

11.19. Outline of proof of Theorem I11.7L We outline the proof from [AN] . Let 

(114) /^P(n)(^) = - E ll^'a^^iriogllP^^lP, 

|a|=?i 

and let 

(115) /,^(„,(^) = - J2 ||PaV^||'log||Pa^||^ 

|a|=ri 

(1) Let K = By the main estimate and the entropy uncertainty inequality, one has 

hrir.,{A) + V")(^/.) > 2(rf - l)n + il^l±^^l^n^(;^) + 0(1). 

(1 - ^) 

(2) Fix no,nE{h). Write n = qrio + r,r < uq. Then by sub-additivity, 

n ~ Uq nEih) riQ 
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(3) As /j, — s> 0, this gives 

WM + kl^^ > 2id - l)n + id-l + 26)X^..^ ^ _ R{no,h) ^ 
no no (l-c") 2no 

(4) Now take the hmit ^ — 0. The expressions | I-Pq^^aI P and the Shannon entropies tend 
to their classical values. Hence, hp(„g) (iph) tends to the KS entropy of the quantum 
limit. Thus, the lower bound implies 

^noil^) ^ ^ _ {d-l + 25)Amax 



no - 2(1-5) 

Since 6, 5' are arbitrary one can set them equal to zero. Then let no — > oo to obtain 
the KS entropy of the quantum limit. 

11.20. Sketch of proof of Theorem 111.51 [X]. The proof of Theorem 111.51 is less struc- 
tural than the proof of Theorem 111.71 There is no magic weapon like the entropic uncertainty 
principle. In its place there is the construction of special covers by cylinder sets and counting 
arguments for the number of elements in the cover. These covering arguments probably have 
further applications and (although quantum) seem more geometric than the argument based 
on the entropic uncertainty principle. 

We mainly discuss the statement that /itop(supp Uq) > ^. We recall ( §11. 6p that an in- 
variant set F is with htop{F) < ^{1 — 6) can be covered by at Ce^^'^'-^"'^-* cylinder sets of 
length A^. We further recall that, by the main estimate f llOSp . each of the cylinder sets has 
quantum measure |/i;i(C)| < (7e-'^/2iogft-AfA/2 uniformly for N < d \ \ogh\ for any We 
want to combine these estimates to show that i^oiF) < 1 for any quantum limit z/q. But we 
cannot simply multiply the measure estimate for C in the cover times the number of cylinder 
sets in the cover. The main estimate being fixed, the only things left to work with are the 
covers. 

As will be seen below, one needs to introduce special covers adapted to the eigenf unctions. 
The covering number of such covers has a sub-multiplicative property. We motivate it by 
going through a heuristic proof from [Aj. We also add a new step that helps explain why 
special covers are needed. See also §11.211 where the pitfalls are explored when F = 7 (a 
closed hyperbolic orbit). 



11.20.1. Covers and times, I. Let F C S be a shift-invariant subset with htop{F) < -^{1 — 5). 
Let Wn C S„ be a cover of minimal cardinality of F by n cylinder-sets. Now introduce a 
new time N » n and define 



(116) ^N{Wn,T) 



{N — cylinders [ao, • • • , o^n-i] such that 



#{je[0,N~n]:[aj,...,aj+„^i]&Wn} y 

N-n+1 — ■ 



They correspond to orbits that spend a proportion > r of their time in Wn- In terms of the 
(ergodic) Liouville measure (transported to S), for almost every point of S*M (or its image 
under /), the proportion of the total time its orbit spends in Wn is the relative measure of 
(the union of the cylinder sets in) Wn- If this measure is smaller than r, then H^iWnjT) 
is a rare event as A^ ^ 00. Its cylinders satisfy cr-^(C) fl Wn = C for a proportion r of 
< J < A^. If one imagines partitioning the elements of Wn into A^ cylinders, then only a 
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(relatively) sparse subset will go into SAr(lV„,r). The next Lemma gives an upper bound 
on that number: 

Lemma 11.13. (see Lemma 2.3.1 of \K\) For all uq there exist n > tiq and Nq so that for 
N > Nq and r G [0, 1], the cover Wn of F of minimal cardinality satisfies, 

where i is the number of elements in the partition. 

If r is chosen sufficiently close to 1, this implies: 

Corollary 11.14. (See (2.3.1) of If F is an invariant set with htop{F) < ^{1-6), 
thcTi 

We now re-write this estimate in terms of the semi-classical parameter h. If r is sufficiently 
close to 1 and e > is sufficiently small, there exists ^9 > so that for N > i)\ log^|, 

(117) ^N(Ml-5))^(l-r)Nil+n) log£ < _ ^) /,<i/2e^A/2 _ 

Indeed, let a = —S^ + (1 — r)(l + n) \ogi. For r very close to 1 it is negative. We then 
choose d so that h~°"^ < (1 - e)h''-/'^, i.e. so that lad] > d/2. 
By the main estimate (11081) and (11171) . we obtain the following 

Corollary 11.15. (see (2.0.1)) 

Forr sufficiently close to 1 and e close to zero, there exists > so that for N > ^\ logh\, 

(118) \fxni^NiWn,T))\<l-e. 
Here, /i/,(S^(W„, r)) means /ih(UceSjv(i^n,r)) '^)- 

11.20.2. Heuristic Proof. We now present a heuristic proof of Theorem 1 1 1 . 5 1 from Section 2 
of ^ that explains the need for some rather technical machinery introduced below. 

We recall that fj,n. is a shift ((T)-invariant measure on S. Let us temporarily pretend that 
/x/j is also a positive measure. Then we would have ior N > {} \ log^|, 

— ( 1 X^N-n-l ^ \ 

< {l-r){l-e) + T. 

Here we used (IllSp . and that the sum of the fit measures of T,N(Wn, t) and its complement 
add to one. Since Wn is fixed, the weak convergence fin. /Uo implies, 

(120) WiWn)\ < (l-r)(l-£)) + r < 1. 

Since F C Wn, the same estimate applies to F. So F cannot be the support of fiQ. "QED" 
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In the above steps, we use that a ^Wn C Siv(W„, r) for A; < — n — 1 and that 

iV-n-l ^ JV-n-l 

(121) Yl 1-'=H.„<1, E K-'=w„<r on J:^{Wn,rr. 

k=0 k=0 

Remark: There are two very innovative ideas in this argument. It is somewhat reminiscent 
of the averaging argument of Theorem 14.11 but there is no averaging over eigenvalues and no 
appeal to ergodicity (since there is no reason why a weak* limit of fin, should be an ergodic 
measure). The first idea is to study the special observables Pa- They are not expressed in the 
form Opn.{a), i.e. in terms of a symbol, and it would be counter-productive to do so. Thus, 
the argument of (11301) takes place entirely on the quantum level and does not make use of 
any properties of the weak* limit (which are unknown); the semi-classical limit is taken after 
the inequality is proved. The key innovation is the third step, where the sets T,N(Wn, r) are 
introduced and the sum is broken up into one over T,N(Wn, t) and one over its complement. 
As is evident from (11211) . EAr(PV„, r)^ arises naturally as an approximation to the r-sublevel 
set of ^k=o^~^ lo-'=VK„- But note again that the the decomposition occurs in Sat and is 
used the on the quantum level. This is very different from using symbols of the C and their 
time averages to to define the corresponding subsets in S*M\ 

The flaw in the argument is in the third line, where we pretended that fin was a positive 
measure. We recall that fin,{C) converges to /io(C), which is a positive measure. But this is 
for a fixed C. But it is also close to being a probability measure for cylinder sets of length 
< k\ logh\. But (11181) is only valid for N > i)\ logh\. Since i) > k, one cannot use them both 
simultaneously (see §11.101 for K,i)) 

One can attempt to get around this problem using the positivity of the underlying state 
Ph- Define the function 

N-n-l 



: = 

N-n 

k=0 



Then in the third line we rigorously have, 

(122) MW:')\ < Pn {iW!:r{W:!))'" = WW^^OpnHl- 



Here, is the quantization as a cylinder set operator. We further denote by T.N{Wn,T) 
the operator equal to the sum of C for C G J2]\r(Wn, t). Then 

So (with Opnix) the energy cutoff), 

(123) \\W^Opn{x)M < \\W^^,AVK,r)OpMM + \\W^^rAwZryOpn{x)M- 
By dm]), 

' \\Wj;'i:^(W^,T)Opn{x)H\ < \\^MW^,r)Opn{x)H\ 
\\Wj:'j:N(wZrYOpn{x)H\ < r\\^NiW^TyOpn{x)H\- 
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To complete the proof of the final equality, need an estimate of the form, 

(124) \\E^{VK,T)Opn{x)H\ <l-0- 
We also need 

(125) \\J:MwZryOpnM < 1 - \\i:,AW^,T)Opn^n\\ + Oih). 

The latter holds formally since SAr(PV„, r), Sjv(PV„, r)'= are semi-classically complementary 
projections. But the time N > -^l logh\ may be too large for this semi-classical estimate. 

The heuristic arguments are sufficiently convincing to motivate a difficult technical detour 
in which a certain sub-multiplicativity theorem is used to reduce the in the argument 
from N > -t^l logh\ (above the Ehrenfest time) to < k| \ogh\ (below the Ehrenfest time). 
It is the analogue in [X] of the sub-additivity step in [ANJ , and its purpose is to allow the 
argument above to use only cylinder sets of length K|log/i| (or balls of radius >> h}^"^) 
where semi-classical analysis is possible. But one has to make sure that the covering number 
estimates in Lemma 111.131 and flll7p do not break down. One has to choose special covers 
so that there are not too many elements in the cover of F. This step is another crucial 
innovation in [X], and in some ways the subtlest one. 

11.20.3. {h, 1 — 0, n)-covers. The sub-multiplicativity step requires a new concept and a new 
parameter 6. Fix {h,n,6) and consider a subset W C 

Definition: Say that W is an {h, 1 — 6,n) — cover of S if 

(126) II CnOpn{x)H\<0- 

Here, is the complement of W\ W is an {h, 1 — 9,n) cover of S if in this quantum sense 

the measure of is < 1 (it is called "small" in [Aj but in the end small means < 1; See the 
Remark at the end of §11.20.41 ). In other words, the union of the cylinder sets of W might 
not cover all of S (think S*M), but the mass of the eigenfunctions ipn. outside of the covered 
part is rather small. In this quantum sense, it is intuitively a union of cylinder sets whose 
measure is > 1 — ^. By definition, if W is an {h, l — 6,n)— cover, then it is an {h, 1 — 6', n) — 
cover for any 6' > 6. 

The covering number of such a cover is defined by 

(127) Nn{n, 9) = : W is a {h,l - e,n) - cover of S}. 
The sub-multiplicative property is 

Lemma 11.16. (sub-multiplicativity) 

(128) Nn {kn, ke{l + 0(nr))) < Nn{n, Of. 

We will not attempt to explain this inequality but rather will point out its consequences. 
One is a lower bound on growth of counting numbers of minimal (/z, 1 — ^, A^)-covers as A^ 
grows (Lemma 2.2.6 of [X]): 

Lemma 11.17. Given 5 > 0, there exists -d so that: 

• (i) For N = 'd\\ogh\, Nn{N,e) > (1 - ^)e^^(i-il); 

• For N = k\ \ogh\, K<^, Nn{N, f 0) > (1 - 6)1 c^^'^^-ts) ; 
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Estimate (ii) follows from (i) and sub-multiplicativity. The time in (ii) is now below the 
Ehrenfest time. The sub-multiplicative property then has the following consequence: 

Lemma 11.18. (See [X], (2.3.2); compare Iil26\) ) Let have the meanings in Hll.lOl Let 

Wn <ZTjn he a cover of F by n-cylinders of minimal cardinality; and let N = k,\ log^|. Then 
for any 6 < 1, 

II 4oMx)^,||>^^. 

CgSjv(VK„,t)= 



Proof. (Sketch) The inequality goes in the opposite direction from that in the definition of 
{h, 1 — 9, iV)-covers in flll.lQp . Hence the Lemma is equivalent to saying that T,j^(W, r) is not 
a (^, 1 — ^6, N)-coveT. This follows by comparing the upper bound in Lemma [H . 1 31 with the 
lower bound in llL17[ If Sjv(W„, r) were such a cover, one would have || log(l — ^^)| > N^, 
which can only occur for finitely many if 6' < 1. The cardinality of J2n{W, t) is too small 
to be {h, 1 — 6, N)-coveT when N = k\ logh\. □ 

The lower bound is a key step in the proof of Theorem 111.51 Intuitively, it gives a lower 
bound on eigenfunction mass outside the set of orbits which spend a proportion r of their 
time in Wn- The lower bound on the mass outside T,N{Wn,T) eventually prohibits the 
concentration of the quantum limit measure on F. 

Using that CfiOpnix) ^^re roughly orthogonal projectors with orthogonal images for distinct 
cylinders C, this leads to: 

Lemma 11.19. (see [X], (2.4.5)) Continue with the notation of Lemma \11.18[ For N = 

k\ log h\, we have 

Eces. \fin{C)\ = l + Oih'), 

The second inequality of Lemma [11.191 again goes in the opposite direction to the defining 
inequality fll26p of a. {h,l — 9, n) cover. 

11.20.4. Completion of Proof. We now sketch the proof of the Theorem (see [X], (2.5.1) - 
(2.5.5)). It is similar to (11191) but circumvents the positivity problem by decreasing N to 
below the Ehrenfest time. 

Proof. (Outline, quoted almost verbatim from \K\ ) 

As above, the quantum measure is shift invariant (see [X], Proposition 1.3.1 (ii).) 

We now fix n and let N = n\ log^|. The crucial point is that this is less than Ehrenfest 
time, unlike the choice of in the first attempt. As above, Wn C S„ denotes an {h, 1 — 6, n)- 
cover of minimal cardinality of F by n-cylinders. In the following, we use the notation 

(129) ^l^,{Wn)■.= Y. ^^^^{C), ^^n{T.N{Wn,T)) := ^in{C). 

c<^w„ ceSjv(w^„,T) 

We then run through the steps of (11191) again; 
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1 sr^N-n-1 
N-n 2^k=0 



( 1 Y^A"-' 



(130) 



< 



< 



< (i-T){i-i^ef)+T + om. 

Above, Lemma [11.191 is used to circumvent the positivity problem in (11191) . We also used 
(again) f ll2ip and the main estimate (11181) to bound the A^-cylinder sets in EAr(W„,r) with 
N = k\ log^l, but with the time n fixed. 

Now let h ^ 0. Since Wn is fixed, the weak convergence fin fiQ implies, 

(131) MWn) < (l-r)(l-(f^^)2)+r. 

Since F C Wn, the same estimate applies to F. Since it holds for all ^ < 1, it follows that 



(132) /^o(F)<(l-r)(l-(^)^)+r = l-,^. 

This completes the proof that supp ^ F . QED 

To obtain the measure statement of Theorem 111.51 let 



r) < 1. 



(133) 



{a;G5*M:/.,(/iS)<^(l-5)}. 



It follows from the Shannon-McMillan theorem that for all a > 0, there exists C I 5 with 
{Xoih\Is) ^ such that can be covered by e"(2"(i~'5+") n-cylinders for sufficiently large n. 
It follows from ffT32l) that 



(134) 

Let a 
(135) 



/io(/r) < (1 
and one has, 



r(5-«)) l-( 



K 



a] 



+ t{6 



a) 



K 



MS*M\Is)>{l-r{6)){j^) 



□ 



11.21. F = 7. To help digest the proof, let us run through the proof in the simplest case 
where F is a single closed hyperbolic geodesic 7 and M is a compact hyperbolic manifold of 
constant curvature —1. We would like to show that z/q 7^ fi^, or equivalently that 1^0(7) < 1- 
What we need to show is that, if z/o(-F) = 1, then it requires a lot of cylinder sets to cover F. 
Hence we try to get a contradiction from the fact that only a few cylinder sets are needed 
to cover 7. In fact, for each n, 7 may be covered by one n- cylinder 

^nil) = {["0, . . . ,ao, . . . ,an]}, 

namely the one specified by the indices of the cells Pa that 7 passes through starting from 
some fixed xq G 7. To obtain a tube around the orbit, we take the union PVn(7) : = 
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Uj=o ''^nil)- We could use the two-sided symbols, going symmetrically forward and 
backward in time to obtain a transverse ball of radius e~"/^. 

First, let us try a naive heuristic argument. The main estimate (11081 ) gives 

(136) \fih{Wn{i))\ < Ch-'^'^e~'"^'^{l + 0(e))", uniformly for n<K\ log/i|, 

for any > 0. It is now tempting to put = (1 + e)|log/i| with e > to obtain 
\iJ^hiWn{h,e))\ ^ But we cannot take the semi-classical limit of this estimate. If we 
did so formally (not rigorously), it would seem to imply that 2/9(7) = 0^ which is stronger 
than the rigorous result 2/9(7) < 1 possibly false in some examples (if weak scarring in 
fact occurs). 

The problem is that the classical set Wn(^fi,e) is a kind of /i-dependent Bowen ball or tube 
around 7 of radius e~^"(^'^) = h}/'^^^ . The uncertainty principle prohibits microlocalization to 
tubes of radius < and (as discussed in it is difficult even to work with neighborhoods 
smaller than h^~^ or at best h^\\ogh\. In particular, we do do not have control over the 
weak limits of \iih{Wn(n,e))\, and there is no obvious reason why it should have semi-classical 
asymptotics. It is not a positive measure on this length scale, so we do not know if it implies 
^^0(7) = 0. On the other hand, if we decrease n below the Ehrenfest time, so that the matrix 
element has a classical limit, then the bound is trivial. 

The strategy of [A] is to use the sub-multiplicative estimate for {h, 1 — 6, A^)-covers of 7 of 
minimal cardinality to decrease the length N below the Ehrenfest time (see Lemma [11.17p . 
We first consider the cover Wn '■= Wn{'y), which is a minimal cover of 7 by cylinders of 
length n. We then form T,N(Wn,T). At least intuitively, Wn corresponds to a Bowen ball of 
radius e""''^, i.e. a transverse cube of this radius in the stable and unstable directions. The 
cylinders of E7v(VF'n, t) satisfy cr-' (C) fl Wn = C for a proportion t of < j < N. This carves 
out a very porous set in the stable/unstable cross section, suggesting why the lower bound 
in Lemma [11.191 should be true. 

It would be interesting to connect the argument of Theorem 1 11. 5 1 in this case to the problem 
in ^ (see §5.ip of estimating the eigenfunction mass inside a shrinking tubular neighborhood 
of 7, i.e. to give an upper bound for the mass {Opfl{xi)'^^J.^ V^/^)) shrinking tube around 7 

or (better) the mass {Opfi{x2)'ffj.y V^/^) ^ shrinking phase space tube. However, the methods 
of [X] apply equally to the cat map setting and in that case, there are eigenfunctions whose 
mass in shrinking tubes around a hyperbolic fixed point have the same order of magnitude 
as in the integrable case up to a factor of 1/2 [FNB] . So unless one can use a property of 
A-eigenfunctions which is not shared by quantum cat map eigenfunctions, one cannot expect 
to prove a sharper upper bound for eigenfunction mass in shrinking tubes around hyperbolic 
orbits in the Anosov case than one has in the integrable case. 

11.22. Some ideas on the main estimate. The main estimate is independently useful. 
For instance, it has been used to obtain bounds on decay rates on quantum scattering [NZj . 
In this section, we sketch some ideas of the proof. We do not go so far as to explain a key 
point, that the estimate is valid for times t < k\ logh\ for any k. 

The goal is to estimate the norm of the operators P*, t/" Pa Op(x*-"'') in terms of h and n. 
Since U{n) is unitary it suffices to estimate the norm of the cylinder set operators 

Kn:=Pa„UP^„^^---UPa,Opix). 
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Let Kn{z,w) be the Schwartz kernel of this operator. By the Schur inequahty, the operator 
norm is bounded above by 



The inequahty only applies when Kn{z,w) G L^{M), but by harmless smoothing approxi- 
mation we may assume this. 

We may get an intuitive idea of the exponential decay of the norm of the cylinder set 
operators as follows. The wave group (or the semi-classical Schrodinger propagator) f/* on a 
manifold without conjugate points can be lifted to the universal cover (see fl25|) ) where it has a 
global Hadamard parametrix f|T71) . This parametrix shows that the wave kernel on the univer- 
sal cover is closely related to the spherical means operator Ltf{x) = /^.^^ f{g^{x,uj)dfix{uj), 
where dfi^ is the Euclidean surface measure on S*M induced by the metric. Both operators 
have the same wave front relation, and indeed the wave front set of U{t,-,y) is precisely 
the set of (co-)normal directions to the distance sphere St{y) centered at y. The operator 
relating [/* and Lt is a dth order (pseudo-differential) operator in {t,x), and in the Fourier 
transformed picture in [2] it gives rise to the factor h~'^^'^ . As time evolves, the surface vol- 
ume of the distance sphere in the universal cover grows exponentially. Since [/* is unitary, 
its amplitude must decay exponentially. In fact this is evident from f[T7l) since 6(x,?/) grows 
exponentially in the distance from x to so 0~^/^(a;, y) decays exponentially. 

The cylinder set operators Cn = Pa„{n) ■ ■ ■ Pa^{l)Pag{0) are given by a sequence of al- 
ternating applications of the propagator U and the cutting off operators Pa (which are 
multiplications). At least heuristically, in estimating the Schur integral for Kn{x,y), one is 
first cutting off to a piece of the distance of the distance sphere Si{y) of radius one centered 
at y. Then one applies U and this piece expands to 5*2 (y). Then one cuts it off with Pa^- As 
this picture evolves, one ends up after n iterates with a piece of the distance sphere Sn{y) 
of radius n centered at y which is has the same size as the initial piece. But the spherical 
means and wave group are normalized by dividing by the volume, so the contribution to the 
Schur integral is exponentially decaying at the rate e~". 

11.23. Clarifications to [X]. 

(1) The 6 on p. 447 is a new parameter < 1 which is implicitly defined above (2.0.1). 

(2) The 6 in (2.3.2) is also implicitly defined by this inequality. Wn-^ is not assumed to 
be a cover. 

(3) In (2.5.1) - (2.5.5), the S in (T~'=S(iy„J should be removed. 

(4) As mentioned above, the notation fih{T,iy(Wn,T)) is not quite defined but means 

/^^(Uc6EAr(iy„,r)) 



This survey has been devoted to the asymptotics of matrix elements {A(fj,ipj). In the 
final section, we address the question, "what applications does the study of these matrix 
elements have?" 

In automorphic forms, the matrix elements are related to L-functions, Rankin- Selberg zeta 
functions, etc. So there is ample motivation for their study in the arithmetic case. Such 
applications lie outside the scope of this survey, so we refer the reader to Sarnak's expository 
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12. Applications to nodal hypersurfaces of eigenfunctions 



RECENT DEVELOPMENTS IN MATHEMATICAL QUANTUM CHAOS 



67 



articles |Sar2l ISar3] for background on and references to the large literature in arithmetic 
quantum chaos. 

Here we are concerned with applications to the classical and geometric analysis of eigen- 
functions on general Riemannian manifolds. We present an application of quantum ergodicity 
to the study of nodal (zero) sets of eigenfunctions on real analytic (M, g) with ergodic geo- 
desic flow. It is not implausible that one could use related methods to count critical points 
of eigenfunctions in the ergodic case. 

The nodal hypersurface of (pj is the zero set jVj = {x G M : (fj{x) = 0}. We are interested 
in its distribution as Xj oo. We study the distribution through the integrals 

[ fdS 

of / over the nodal line with respect to the natural surface measure. If / = Ijy is the 
characteristic (indicator) function of an open set U C M, then 

(137) / luds = W-^{U nUj) 

is (^ — 1) dimensional Hausdorff measure of the part of Mj which intersects U . When U = M 
the integral fll37p gives the total surface volume of the nodal set. The principal result on 
volumes is due to Donnelly- Fefferman |DFj . proving a conjecture of S. T. Yau |Ylj : 

Theorem 12.1. |DFj (see also |Linj ) Let {M,g) be a compact real analytic Riemannian 
manifold, with or without boundary. Then there exist Ci,C2 depending only on {M,g) such 
that 

(138) ciA < n^-\Z^J < C2X, {A<f^ = AVa;ci,C2 > 0). 

The question is whether one can obtain asymptotics results by imposing a global dynamical 
condition on the geodesic flow. 

The hypothesis of real analyticity is crucial: The analysis is based on the analytic contin- 
uation of eigenfunctions to the complexification Mq of M. The articles |DF^ [Lm J also used 
analytic continuation to study volumes of nodal hyp ersurf aces. It is the role of ergodicity in 
the complex setting that gives strong equidistribution results. 

12.1. Complexification of {M,g) and Grauert tubes. By a theorem of Bruhat and 
Whitney, a real analytic manifold M always possesses a unique complexification Mc , an open 
complex manifold in which M embeds l : M ^ Mc as a totally real submanifold. Examples 
are real algebraic subvarieties of which complexify to complex algebraic subvarieties 
of C^. Other simple examples include the complexification of a quotient such as M = 
M™'/Z''" to Mc = C™/Z''". The complexification of hyperbolic space can be defined using the 
hyperboloid model H" = {xj + ■ ■ ■ x"^ — xl^^^ = — 1, x„ > 0}. Then, = {{zi, . . . , ^n+i) £ 

The Riemannian metric determines a special kind of plurisubharmonic exhaustion function 
on Mc |GSll IGS2t ILemSlj . namely the Grauert tube radius ^/p = ^/p^ on Mc, defined as 
the unique solution of the Monge- Ampere equation 



{ddy/p)"" = SMt^dVg, t*{iddp) = g. 
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Here, SM^,dVg is the delta-function on the real M with respect to the volume form dVg, i.e. 

/ — * Juf'^^g- shown in |GSlj that y/p{C) = iyi^XTC) where r^{x,y) is the squared 
distance function in a neighborhood of the diagonal in M x M. 
The Grauert tubes Mr is defined by 

Mr = {CeMc: VpiO < r}- 

The maximal value of tq for which y^p is well defined is known as the Grauert tube radius. 
For T < To, Mr is a strictly pseudo-convex domain in Mq. Using the complexified exponential 
map (x,^) — >■ expxi^ one may identify Mr with the co-ball bundle i?*M to Mr. Under this 
identification, ^/p corresponds to \^\g. The one-complex dimensional null foliation of dd^/p, 
known as the 'Monge- Ampere' or Riemann foliation, are the complex curves t + ir T'y{t), 
where 7 is a geodesic. The geometric Grauert tube radius is the maximal radius for which 
the exponential map has a holomorphic extension and defines a diffeomorphism. 

12.2. Analytic Continuation of eigenfunctions. The eigenfunctions are real analytic 
and therefore possess analytic continuations to some Grauert tube M^ independent of the 
eigenvalue. To study analytic continuation of eigenfunctions, it is useful to relate them to 
analytic continuation of the wave group at imaginary time, U{iT,x,y) = e~'^^{x,y). This 
is the Poisson operator of {M,g). We denote its analytic continuation in the x variable to 
by Uclir, (,y). In terms of the eigenfunction expansion, one has 

00 

(139) U{tT,C,y) = Y,e~^'^iff{Ov,{y), {C,y)eM,xM. 

j=0 

To the author's knowledge, the largest e has not been determined at this time. In partic- 
ular, it is not known if the radius of analytic continuation of the wave kernel and of the 
eigenfunctions is the same as the geometrically defined Grauert tube radius. 
Since 

(140) Uc{it)vx = c-^Va, 

the analytic continuability of the Poisson operator to Mr implies that every eigenfunction 
analytically continues to the same Grauert tube. It follows that the analytic continuation 
operator to Mr is given by 

(141) Ac{T) = Uc{tT)oe^^. 

Thus, a function / G C°°(M) has a holomorphic extension to the closed tube ^/p{C) < r if 

and only if / G Dom{e^^), where is the backwards 'heat operator' generated by \fK 
(rather than A). 

Let us consider examples of holomorphic continuations of eigenfunctions: 

• On the flat torus M^/Z™, the real eigenfunctions are cos(fc, x), sin(A;, x) with k G 
2t^7/^. The complexified torus is <C^/'L'^ and the complexified eigenfunctions are 
cos(/c, C), sin(fc, Q with = x + i^. 

• On the unit sphere S*™, eigenfunctions are restrictions of homogeneous harmonic 
functions on The latter extend holomorphically to holomorphic harmonic 
polynomials on C"^"^^ and restrict to holomorphic function on S"^. 
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• On H™, one may use the hyperbolic plane waves e^*'^^^)^^'''^, where (z, b) is the (signed) 
hyperbolic distance of the horocycle passing through z and b to 0. They may be 
holomorphically extended to the maximal tube of radius tt/A. See e.g. |BRj . 

• On compact hyperbolic quotients H"^/r, eigenfunctions can be then represented by 
Helgason's generalized Poisson integral formula, 

Mz) = [ e(^^+i)<^'^>rfr,(6). 
Jb 

To analytically continue (fx it suffices to analytically continue {z, b). Writing the latter 
as (C, "we have: 

(142) y.^(C) = / e(*^+i)«'''>drA(6). 

JB 

• In the case of ID Schrodinger operators, several authors have studied analytic con- 
tinuations of eigenfunctions and their complex zeros, see e.g. |Hez^ lEGSj . 

12.3. Maximal plurisubharmonic functions and growth of tp^. In the case of domains 
VL C C"^, the maximal PSH (pluri-subharmonic) function (or pluri-complex Green's function) 
relative to a subset E dVLi?, defined by 

Ve{C) = snp{u{z) : u G PSH{Q),u\e < 0,u\an < 1}. 

An alternative definition due to Siciak takes the supremum only with respect to polynomials 
p. We denote by the space of all complex analytic polynomials of degree N and put 
= {p G : Iblli, < 1, I bib < e}. Then define 

log<l>f(C) =sup{4logb7v(C)l -peV^}, log^ij = limsuplog$f. 

Here, ||/||k = sup^g^ l/(-2)|- Siciak proved that log^E = Ve- 

On a real analytic Riemannian manifold, the natural analogue of is the space 

'H^ = {p= ^ ajV^A,, ai, . . . ,aAr(A) G M} 

j:\j<X 

spanned by eigenfunctions with frequencies < A. Rather than using the sup norm, it is 
convenient to work with based norms than sup norms, and so we define 

7V(A) 

i:Aj<A j=l 

We define the A-Siciak approximate extremal function by 

^Uz) = sup{|^(z)|VA: ^ e Hx; Mm ^ 1}, 
and the extremal function by 



<I>m{z) = snp<^lj{z) 

A 



It is not hard to show that 



(143) $m(^) = lim \ logn[o,A(C, = Vp- 

A— >oo A 
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12.4. Nodal hypersurfaces in the case of ergodic geodesic flow. The complex nodal 
hypersurface of an eigenfunction is defined by 

(144) Z^c = {C e BIM : (^^(C) = 0}. 

There exists a natural current of integration over the nodal hypersurface in any ball bundle 
B^M with e < eo , given by 

(145) {[Zr],^) = ^f dd\og\^l\^A^= f ^eD(™-i'-i)(5:M). 

^ ^TT Jb*M Jz c 

In the second equality we used the Poincare-Lelong formula. The notation T>'^'^~^''^~^\B*M) 
stands for smooth test (m — 1, m — l)-forms with support in B^M. 

The nodal hypersurface Z^c also carries a natural volume form \Z^c\ as a complex hyper- 

^m — 1 

surface in a Kahler manifold. By Wirtinger's formula, it equals the restriction of , to 



m — 1 

(m-1)! 



Theorem 12.2. Let (M, g) he real analytic, and let {(pj^} denote a quantum ergodic sequence 
of eigenf unctions of its Laplacian A. Let {B*^^M, J) be the maximal Grauert tube around M 
with complex structure Jg adapted to g. Let e < cq. Then: 

-^[Z r ] ^ -dd^ weakly in V'^'^'^\B*M), 
in the sense that, for any continuous test form e 'D^"^~^'"^~^^{B*M), we have 

^Jk Jz c ^ JBiM 



Equivalently, for any e C{B*M), 



KJzc ("^-1)! t^Jb*m (m-1)! 



Corollary 12.3. Let {M,g) be a real analytic with ergodic geodesic flow. Let {(fiji^} denote 
a full density ergodic sequence. Then for all e < eg, 

■^[Z^c ] -dd^, weakly in V'^^^^\B*M). 

Remark: The limit form is singular with respect to the Kahlerform iddp on Mq adapted to 
the metric. Its highest exterior power is the delta function on the zero section (i.e. the real 
domain M). 

In outline the steps in the proof are; 

(1) By the Poincare-Lelong formula, [Z^c] — iddlog\(p'^\. This reduces the theorem to 
determining the limit of j^og \ (p'x\- 

(2) Iv^^P, when properly L^ normalized on each dMr is a quantum ergodic sequence on 
dMr. 
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(3) The ergodic property of complexified eigenf unctions implies that the norm of |<y5^p 
on dM-r is asymtotically y^. Thus, normalized log moduli of ergodic eigenfunctions 
are asympotically maximal PSH functions. They are A-Siciak extremal functions. 

(4) i log qp'^''^ > for an ergodic sequence of eigenfunctions. Hence, 

(5) {dd\og\^'^\ idd^, 
concluding the proof. 

Remark: We expect a similar result for nodal lines of Dirichlet or Neumann eigenfunctions 
of a piecewise smooth domain f2 C with ergodic billiards. We also expect similar results 
for nodal lines of such domains which touch the boundary (like the angular 'spokes' of 
eigenfunctions of the disc). In |TZ1] . we proved that the number of such nodal lines is 
bounded above by CqX. This is rather like the Donnelly-Fefferman upper bound. However 
the example of the disc shows that there is no lower bound above zero. That is, there 
are sequences of eigenfunctions (with low angular momentum) with a uniformly bounded 
number of nodal lines touching the boundary. We believe that ergodicity is a mechanism 
producing lower bounds on numbers of zeros and critical points. It causes eigenfunctions to 
oscillate uniformly in all directions and hence to produce a lot of zeros and critical points. 
We expect this to be one of the future roles for ergodicity of eigenfunctions. 

Remark: Precursors to Theorem 112. 21 giving limit distributions of zero sets for holomorphic 
eigen-sections of Kahlerquantizations of ergodic symplectic maps on Kahlermanifolds can be 
found in [NV2^ ISZj . In that setting, the zeros become uniformly distributed with respect to 
the Kahler form. In the Riemannian setting, they become uniformly distributed with respect 
to a singular (1, 1) current relative to the Kahler form ddp on Mc- 
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